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1.0  INTRODUCTION 


This  report  is  the  final  report  on  the  research  program  on  the  fundamental  processes 
in  partially  ionized  plasmas  conducted  in  the  High  Temperature  Gasdynamics  Laboratory  at 
Stanford  University.  This  research  is  supponed  by  a  grant  from  the  Air  Force  Office  of 
Scientific  Research  (AFOSR-88-0264)  and  is  currently  conducted  under  the  direction  of 
Professor  Charles  H.  Kruger.  Five  graduate  students  have  completed  their  Ph.D.  under 
this  program.  Four  of  these  currently  have  faculty  positions  at  other  universities. 

Several  space  power  and  propulsion  systems  of  potential  long  range  interest  to  the 
Air  Force  involve  partially  ionized  plasmas.  Such  systems  include  MPD  thrusters,  both 
open  and  closed  cycle  MHD  power  generation,  high  velocity  reentry  and  thermionic  energy 
conversion.  Although  the  specific  configurations,  the  exact  operating  conditions,  and 
which  of  the  competing  systems  will  prove  to  be  most  useful  in  the  long  term  remain  to  be 
established,  it  is  imponant  at  this  time  to  provide  a  broad  fundamental  research  base  in 
support  of  development  activity.  In  particular,  there  are  a  number  of  key  issues  regarding 
the  properties  and  discharge  behavior  of  partially  ionized  plasmas  and  the  interaction  of 
discharges  with  fluid  dynamics  that  need  to  be  understood  before  the  potential  and 
limitations  of  competing  systems  can  be  fully  evaluated.  In  addition,  it  is  important  that 
outstanding  young  applied  scientists  be  educated  in  these  areas. 

Our  research  on  partially  ionized  plasmas  was  initiated  under  grant  AFOSR-83- 
0108  and  focused  on  three  major  areas: 

1 .  plasma  properties, 

2.  discharge  effects:  plasma  electrode  interaction, 

3 .  interaction  of  discharges  and  fluid  dynamics. 

Extensive  reports  on  areas  (2)  and  (3)  have  been  given  in  previous  Annual  Scientific 
Reports  and  in  the  Final  Report  for  grant  AFOSR-83-0108.  Recent  research,  under  grant 
AFOSR-88-0264,  has  emphasized  area  ( 1 ),  as  does  this  report. 

Results  on  air  plasma  properties  are  described  in  the  following  section. 
Publications  and  presentations  resulting  from  this  work  are  cited  in  Section  3  and  Section  4 
lists  the  personnel  who  have  contributed  to  this  report. 
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:.0  OPTICAL  DIAGNOSTICS  AND  RADIATIVE  EVUSSION  OF  AIR  PLASMAS 


1.  Summary  and  Overview 

At  high  temperatures,  air  is  a  complex  mixture  of  atoms,  moiecuies  and  charged  pardcies: 
commonly  referred  to  as  a  plasma.  The  radiative  properties  of  air  plasmas  are  of  primary  interest 
in  a  wide  range  of  applications,  including  chemical  and  material  processing  (plasma  spray¬ 
coating),  design  of  heat-shields  for  space  vehicles  (atmospheric  re-entry,  aeroassisted  orbital 
transfer  vehicles),  and  detection  of  missile  signatures  (SDI).  Over  the  past  three  decades,  space 
programs  have  stimulated  intense  research  on  air  radiation,  both  experimentally  (see  for  example 
Page  and  Arnold,  1964,  Nerem,  1964,  Schieiber  et  al.,  1973,  Ogurtsova  et  al.,  1968,  Devoto  et 
al..  1978)  and  theoretically  with  the  development  of  numerical  codes  (Kivel  and  Bailey,  1957, 
Meyerott  et  al.,  1960,  Breene  et  al.,  1963,  Allen.  1966.  Whiting  et  al..  1969,  Park.  1985a,  i985b). 
The  projected  manned  missions  to  Mars  have  recently  provided  a  new  thrust  for  radiative  sradies 
as  the  Apollo-shaped  vehicles,  frequently  considered  for  future  atmospheric  manned  reentries,  are 
expected  to  undergo  80%  of  the  stagnation  point  peak  heating  from  shock-layer  radiation  at  the 
design  re-entry  velocity  of  14  km/s  (Tauber  et  al.  1992). 

Air  plasma  modeling  is  a  complex  task  that  requires  the  understanding  of  fundamental 
molecular  processes  to  predict  the  population  of  every  rotational,  vibrational  and  electronic  level 
for  each  species  present  in  the  gas,  along  with  emission  intensities  for  transitions  between  these 
levels.  Modeling  can  be  separated  into  three  modules: 

•  flow  field  module,  to  predict  the  chemical  composition  of  the  plasma. 

•  excitation  module,  to  predict  the  electronic,  vibrational,  and  rotational  level  populations. 

•  radiative  module,  to  predict  the  intensity  emitted  by  ail  possible  radiative  transitions. 

Among  the  various  codes,  NEQAIR  (Park,  1985a  and  1985b)  has  been  widely  used  because 

it  includes  both  an  excitation  and  a  radiation  module,  and  because  these  modules  comprise  most 
physical  processes  of  importance.  Moreover,  NEQAIR  has  been  the  subject  of  ongoing 
development,  including  recent  improvements  in  performance  (speed  increase  achieved  through 
optimization,  vectorization.  and  restructuring  of  the  code)  and  flexibility  (Moreau  et  al..  1992). 

While  experimental  data  are  necessary  to  assess  the  quality  of  the  codes,  they  are  often 
difficuit  and  expensive  to  obtain.  Nevertheless,  in  addition  to  the  previously  mentioned 
experiments  and  to  those  referenced  by  Park  (1985a),  spectral  measurements  have  been  recently 
conducted  in  shock  tubes  (Sharma  et  al.,  1991a.  Sharma  and  Gillespie,  1991b)  and  in  the  bow 
-hock  of  bailisdc  missiles  lErdman  et  al..  1992).  However,  to  numerically  reproduce  these 
experiments,  the  chemical  composition  and  the  various  characteristic  temperatures  (translational. 
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rotanonai.  vibrationai.  and  eiectromc )  of  the  piasma  must  be  obtained  from  computadons.  It  is 
not  clear  whether  the  observed  discrepancies  that  may  arise  between  the  experiments  and  the 
computations  are  then  due  to  errors  in  the  flowtield  computation,  or  in  the  excitation  module,  or 
in  the  radiative  module.  In  order  to  check  the  validity  of  each  module,  a  systematic  approach  is 
therefore  essential  In  particular,  the  radiative  module  must  be  thoroughly  checked  so  that  non¬ 
equilibrium  experiments  enable  investigators  to  check  the  flowfield  and  excitation  modules. 

TIte  work  presented  in  this  repon  has  been  mainly  dedicated  to  providing  a  benchmark 
experiment  to  test  the  radiative  module  included  in  NEQAIR.  and,  based  on  tnese  experimental 
results,  to  improving  the  radiative  model.  For  clarity,  the  original  code  is  termed  NEQAIR- 1,  and 
the  enhanced  version.  NEQAIR-2.  The  adopted  strategy  is  presented  in  the  following  overview. 

In  Sections  2  and  3,  the  plasma  torch  facility  and  the  experimental  device  used  for  the 
emission  measurements  are  desen'oed  in  detail 

In  Section  4,  the  thermodynamic  state  of  the  air  piasma  generated  in  the  torch  is  investigated 
using  emission  measurements  to  determine  eiectromc  and  rotational  temperatures  as  well  as 
electron  number  densities.  The  chemical  composition  is  also  studied  by  means  of  the  kinetic  code 
CHEMKIN  with  an  extended  thermodynamic  data  base  to  generate  reverse  reaction  rates  above 
5,000  K.  The  results  presented  in  this  section  show  that  the  piasma  is  approximately  in  local 
thermodynamic  equilibrium  (LTE)  over  most  of  the  piasma  extent.  An  LTE  temperature  profile  is 
measured  to  an  accuracy  of  a  few  percent  (typically  less  than  3%).  This  single  profile  fully 
characterizes  the  thermodynamic  state  of  the  piasma,  and  provides  the  input  data  to  perform  the 
numerical  simulations  of  Section  5. 

In  Section  5.  the  focus  turns  to  the  testing  of  the  radiative  module  of  NEQAIR  through  a 
comparison  between  the  experimental  emission  spectrum,  measured  and  calibrated  in  intensity 
over  the  spectral  range  2,(XX)  to  8,0(X)  A,  and  numencal  simulations  using  the  NEQAIR  code. 
The  experimental  spectrum,  measured  along  a  plasma  diameter,  includes  contribution  of  regions  at 
temperatures  up  to  7,5(X)  K  which  is  the  temperature  at  the  center  of  the  axisymmetric  plasma 
torch.  This  temperature  range  is  very  well  suited  to  investigate  the  radiation  of  the  major 
radiative  bands  of  air  plasmas.  As  will  be  seen  in  that  section,  the  spectrum  contains  strong 
atomic  and  molecular  features  due  to  NO,  N^,  NT,  CN,  O^,  O,  N,  and  C.  Since  the  plasma 
composition  and  temperature  are  well  known  from  the  previous  section,  it  is  possible  to  simulate 
the  spectrum  numerically  using  the  radiative  module  of  NEQAIR.  The  comparison  with  the 
original  version  of  NEQAIR  revealed  that  some  features  of  the  code  could  be  modified  or 
cnhancea.  Vr.z  revisions,  detailed  in  the  repon.  include  a  thorough  uodate  of  molecular  radiative 
transition  probabilities  and  spectroscopic  constants,  and  the  addition  of  several  NO  band  systems 
(NO  Delta.  Epsilon.  Beta  prime  and  Gamma  prime).  Since  the  C  state  from  which  the  NO  Delta 


Tansidon  oneinates  is  preaissociated.  a  simpiined  coiiision-premsscx:iadon  model  for  mis  state 
has  been  added  to  the  code.  The  effects  of  these  changes  are  discussed,  and  a  final  comparison 
between  the  experimental  spectrum  and  the  results  of  the  ennanced  NEQAIR  code  is  presented. 

Finally,  Secdon  6  presents  measurements  of  the  total  radiadve  source  strength  of  LTE  air 
between  4000  and  7500  K.  Tnese  results  are  then  compared  to  the  predicdons  of  several  codes, 
in  pardcular  the  modified  NEQAIR  code. 
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2.  Experimental  facility:  Plasma  torch. 


High  temperature  plasmas  may  be  conveniently  produced  in  arcs  and  torches.  These  various 
devices  are  described  in  details  by  Dresvin  et  ai  (1977),  Davies  and  Simpson  ( 1979).  and  Boulos  et  al 
(1985).  Plasma  torches  offer  advantages  over  arcs  for  they  do  not  include  electrodes,  hence 
eliminating  contaminadon  problems.  Induction  plasma  torches  may  produce  either  high  or  low 
enthalpy  plasma  flows.  High  enthalpy  torches,  which  inherently  operate  at  low  efficiency,  are  suited 
for  aerospace  applications  such  as  supersonic  re-entry  simulanon  (Deger  et  al.,  1990).  In  contrast, 
low  enthalpy  torches  such  as  the  one  available  in  our  laboratory  are  typically  designed  for  material 
and  chemical  processing  (plasma  spray-coating).  Torches  of  the  latter  type  achieve  efficiencies  on 
die  order  of  50%  and  produce  relatively  uniform  enthalpy  profiles  at  the  nozzle  exit 

Energy  is  coupled  into  the  plasma  by  either  an  oscillating  eiectnc  field  (capacitive  coupling;  or 
an  oscillating  magnetic  field  (inductive  coupling).  The  pnneipies  of  inductive  coupling  are  presented 
in  Section  2.1  since  the  torch  used  in  the  present  work  is  of  that  type.  The  characteristics  of  the 
acmai  torch  are  then  presented  in  Section  2.2.  and  an  estimate  of  the  plasma  properties  in  the 
induction  region  is  given  in  Section  2.3. 


2.1.  Induction  heating 


The  basic  design  for  inductively  coupled  plasma  iICP)  torches  has  not  changed  much  since 
•Jieir  first  introduction  by  Reed  ( 1961).  Gas  is  injected  at  the  bottom  of  a  quanz  tube  surrounded  by 
a  coaxial  induction  coil  traversed  by  a  radio  frequency  currenL  The  rf  current  produces  an  oscillating 
.Lxial  magnenc  field  that  forces  the  free-eiectrons  to  spin  in  a  radial  plane  and  thereby  generates  eddy 
currents.  The  oscillating  magnetic  field  induces  an  oscillating  electric  field.  £(r),  as  may  be  seen 
from  Maxweii's  equation: 


vxe=-M 

dt 

Inteeratins  ever  the  area  shown  in  Fig.  2.1.  and  using  Stokes  theorem. 


(2.1) 


il(VxE).ndA  =  {.E  df  =iX^  =idA. 
Since  the  eiecuic  field  is  tangential  (Gauss  theorem).  (2.2)  reduces  to 


^  nEa  ,  _  nB. 

x(r-h  tir  )i  Ea  -1 - dr  |-  Krt.^  - - -.rar . 

L  dr  j  dr 


and  neglecting  high  order  derivatives. 


(2.2) 


..jj 
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(2.4) 


_  oEoi  aB. 

Eq-rr— =  -r— 
or  at 

This  induced  electromagnetic  force  accelerates  the  free  electrons,  thereoy  increasing  their 
kinetic  energy  e,  or  equivalently  their  transiationai  temperature.  Tg  =  2/3  e/k.  The  tree  electrons  then 
nansier  their  kinetic  energy  to  heavy  particies  through  collisions.  When  ionizing  collisions  produce 
enough  electrons,  the  plasma  becomes  self-sustained. 


Z  ♦ 

I 


I 
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Figure  2.1.  Area  of  integration. 

If  there  were  no  plasma  in  the  coil  region,  the  magnetic  field  radial  profile  would  be  almost 
uniform  (truly  uniform  tor  a  coil  of  infinite  length).  However,  when  the  plasma  is  present,  the  emf 
(hat  acts  on  the  eddy  currents  creates  an  adverse  magnedc  field  (Lenz  law)  that  mcreases  in  intensiu’ 
iowaras  the  center  of  the  torch  and  tenas  to  cancel  the  driving  B  field.  Tne  piasma  therefore  shieids 
electromagnetic  waves  and  limits  their  penetration  to  the  outer  radial  region  of  the  discharge.  This 
phenomenon  was  experimentally  observed  by  Ecken  1 1971)  as  shown  in  Figure  2.2. 

Freeman  and  Chase  ( 1968)  calculated  the  H-field  distribution  in  the  approximation  of  the  so- 
cailed  Channel  model,  anc  obtained  the  following  expression  at  radius  r; 


ber{'j2rld)  +  j  beii-Jlr  I d] 

H(r)  =  2H^ - r7=-l - TP-T- 

her(V2p)  +  jbeiiyPZq) 

where  her  and  bei  are  the  real  and  imaginary  pans  of  the  zero  order  Bessel  funcnon.  Hi  the 
.uTiDiitude  c:  the  magnetic  field,  a-nd  p  a  reduced  parameter.  The  so-called  skin  deoth  parameter  d. 
defined  as 
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(2.6) 


dB.Jl/K\if^CSJ-  \ 

;haracienzes  the  depth  of  penetration  of  the  electromagnetic  wave  into  the  piasma. 


Figure  1.2.  Radial  profiles  of  the  magnetic  field  intensity  at  the  mid-section  of  the  coil  region,  in  the  absence  (a) 
and  presence  (b)  of  the  discharge  (after  Eckert,  1971). 

This  parameter  depends  on  the  frequency  /  of  the  discharge  and  on  the  electrical  conductivity  a  of 
the  gas.  For  an  atmospheric  pressure  air  plasma  at  an  average  temperature  of  8.000  K.  the  electrical 
jonducnviry  is  a=800  ohm'‘rrr‘  (Dresvin  et  ai..  p  42.  1977).  .At  an  operanng  freauency  of  4  MHz, 
the  SKin  aeptn  can  therefore  be  approximatea  by: 

d  “  (4;r*x  10~'  <  800x4  x  10^ )  ^  9  mm .  (2.7) 

which  shows  that  the  electromagnetic  field  penetrates  about  a  quaner  of  the  way  into  our  70  mm 
diameter  plasma  torch.  The  inductive  heating  process  therefore  takes  place  primaniy  in  an  outer  ring 
and  the  central  pan  of  the  piasma  ts  subsequently  heated  by  convection  and  thermal  conduction  of 
heavy  particles.  This  results  in  lower  temperatures  at  the  center  the  discharge  and  in  an  M-shaped 
temperature  profile  in  the  coil  region.  The  skin  depth  could  be  increased  by  decreasing  the  operating 
frequency,  out  this  would  also  result  in  a  decrease  of  the  coupling  efficiency  which  is  deterrmned  by 
•he  couoiing  parameter  k  =  v'2  f.  jd.  where  r„  is  the  torch  radius. 
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-.2.  Description  of  the  piasma  torch 


The  aevice  used  in  this  work  is  a  T.AFA  Model  66  inducnveiy  coupied  piasma  torch  powered 
by  a  radio  frequency  LEPEL  Model  T-50-3  power  supply  operating  at  4  MHz.  This  power  supply 
can  deliver  up  to  120  kV A  or  line  power  to  the  oscillator  plates  with  a  maximum  of  12  kV  DC  and 
-.5  A.  The  oscillator  plates  have  in  turn  a  maximum  RF  power  output  of  50  kW.  The  heat 
discharged  into  the  torch  and  generator  cooling  systerr  is  measured  by  a  combination  of 
thermocouples  at  the  inlet  and  oudet  of  the  generator  and  turbine  flowmeters  in  the  flow  lines.  After 
subtracung  the  heat  lost  into  the  oscillator  tube,  this  information  is  used  to  determine  the  power 
deposited  into  the  plasma. 

The  torch  schematic  (Figure  2.3')  shows  the  5-tum  copper  induction  coil  (mean  diameter 
S.6cm)  encapsulated  between  a  confinement  tube  (inner  diameter  7.6  cm,  thickness  3  mm)  and  a 
Teflon  Dody.  The  confinement  tube  is  made  of  quartz,  a  matenai  that  offers  good  thermal  resistance 
and  is  reiativeiy  transparent  to  eleciromagnenc  fields.  The  outer  Teflon  body  acts  as  an  electrical 
insulator  and  as  an  electrom.''<menc  screen.  The  coil  is  cooled  with  de-ionized  water  to  prevent 
arcing  between  hs  turns,  .although  water-cooling  is  technically  more  complicated  to  operate  than 
air-cooling,  a  higher  power  operating  range  and  thus  higher  heating  of  the  gas  is  achieved. 

The  piasma  torch  can  operate  with  a  variety  of  gases  (argon,  air,  hydrogen,  nitrogen,  methane, 
etc.)  whose  flowrates  are  monitored  in  an  independent  free  standing  panel.  While  the  work  reponed 
here  was  done  with  air,  argon  was  used  during  the  torch  starting  phase,  and  hydrogen  was  added  for 
electron  number  density  measurements. 

The  gas  is  injected  in  axial,  radial  and  swirl  mexies  through  the  manifold  located  at  the  bottom 
'll  the  torch.  .Axial  injection  provides  bulk  movement  to  the  gas  and  is  therefore  useful  when  starting 
'.he  torch.  In  normal  operauon.  only  swirl  and  radial  injectors  are  used,  in  the  rauo  i/3:2/3.  Radial 
mjecuon  provides  the  desired  axisymmetne  pattern  of  gas  flow.  Swirl  injection  creates  a  cold  gas 
sheath  close  to  the  quanz  coniinement  tube  to  prevent  overheating  and  to  reduce  heat  losses  to  the 
cooling  water.  In  addition,  the  swirl  creates  a  low  pressure  zone  in  an  annular  shell  that  favorizes 
gas  recirculadon  in  the  induction  region  and  stabilizes  the  plasma.  For  this  reason,  such  plasmas  are 
called  vonex -stabilized  plasmas.  The  recu'cuiadon  zone  has  been  evidenced  by  several  invesngators. 
expenmentally  by  Chase  (1971)  and  Klubnikin  (1975),  and  theoretically  by  Rhodes  and  Keefer 
( 1989). 

The  water-cooled  hollow  probe  at  the  center  of  the  manifold  is  designed  for  powder  or  gas 
injection  in  commercial  operation.  In  the  present  case,  however  it  is  only  used  when  starting  the 
;orcn  to  carr/  the  free  eieccons  generated  outside  the  torch  body  'oy  a  high  freauency  tesla  ccii 
(Electro-Technic  model  BD-IO-A).  These  free  electrons  are  necessary  to  tngger  the  ionizing 
collisions  that  will  finally  provide  enough  electrons  to  sustain  the  oiasma.  In  this  starting  phase. 
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argon  is  preierrea  over  other  gases  because  us  low  specinc  heat  permits  a  rapid  temperature  increase 
leading  to  ionizadon.  and  because  its  low  thermal  conductivity  prevents  the  hot  zone  trom  cooling  by 
diffusion,  which  would  result  in  thermal  pinching.  Once  the  argon  plasma  is  established,  the  plate 
power  is  increased  and  air  is  progressively  substituted  to  argon. 


The  piasma  then  expands  into  ambient  air  through  a  converging  copper  nozzle.  Different 
nozzle  diameters  from  a  few  millimeters  to  7  cm  can  be  mounted  on  the  torch:  for  the  present  work, 
a  5  cm  nozzle  was  used..  All  measurements  presented  in  this  report  were  conducted  1  cm 
iownstream,  cf  the  nozzle  exit.  In  this  radial  plane,  t.he  moiximum  veiocitx'  is  about  20  m/s  and  the 
maximum  temperature  is  7.500  K.  corresponding  to  a  density  p  of  2.7x10'^  kg.m-^  and  a  dynamic 
viscosity  u  =  i, 60x1  O'*  kg.i  '.s  ’.  With  these  values,  and  based  on  a  nozzle  diameter  of  5  cm.  the 
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Reynolds  number  at  the  measurement  pomt  is  about  i70.  Tne  plasma  let  is  therefore  iaminar  in  the 
measurement  region:  it  remains  so  downstream  of  the  nozzie  over  a  few  centimeters,  then  turns  into 
a  bushy  piume  approximately  50  cm  nigh. 


2  J.  Properties  of  the  air  plasma  in  the  induction  re^on 

The  induction  zone  could  not  be  optically  accessed  here,  and  therefore  no  measurements  were 
made  in  that  region.  However,  an  estimate  of  the  plasma  properties  inside  the  coil  region  will  be  of 
importance,  particularly  in  Section  4.  when  smdying  possible  chemical  non-equilibrium  effects. 
Several  authors  have  investigated  the  induction  region,  both  theoretically  and  experimentally.  We 
present  below  some  of  their  observations  and  how  they  may  oe  applied  to  the  present  torch. 

As  mentioned  in  the  previous  section,  the  center  of  the  piasma  m  the  induction  region  is 
expected  to  be  colder  than  the  annular  shell  into  which  most  of  the  electromagnetic  energy  is 
coupled.  The  temperature  difference  between  the  coupling  shell  and  the  center  depends  on  the  type 
of  gas.  the  skin  depth,  and  the  dimensions  of  the  torch.  This  difference  was  however  reported  to  be 
small,  with  a  maximum  of  1000  K  in  argon  plasmas  (Gol'dfarb  and  Drcsvin.  1965),  and  a  few 
hundred  degrees  for  oxygen,  nitrogen  or  air  plasmas  (Dresvin  et  al.,  p  194,  1972). 


2-9 


Temperature  measurements  were  presented  by  Goykhman  et  ai.  ( 1970)  for  an  air  plasma 
prtxiuced  by  a  torch  simiiar  to  the  present  one;  an  inductive  plasma  torch,  with  a  60  mm  diameter 
quartz  mbe.  a  4-tum  inducdon  coil,  and  50  kW  nominal  power.  The  5.28  MHz  generator  ftequency 
was  slightly  higher  than  in  the  present  torch,  but  Dresvin  et  al.  (1972)  observed  that,  for  argon 
plasmas  at  least,  the  maximum  temperature  in  the  coil  region  is  independent  of  the  operating 
frequency.  Fig.  2.4  shows  a  plot  of  the  temperature  averaged  over  a  plasma  radius  in  the  mid-secdon 
of  the  coil  region  as  a  function  of  the  dissipated  power  and  for  two  flow  rate  cases.  Since  our 
experiments  were  conducted  at  a  flow  rate  of  95  slpm  and  a  deposited  power  of  about  27  kW,  the 
average  temperature  in  the  induction  region  should  not  be  above  8000  K.  and  taking  into  account  the 
observations  from  the  previous  paragraph.  9,000  K  is  a  conservative  estimate  for  the  maximum 
temperature  in  the  induction  region. 
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3.  Experimental  facility:  Optical  measurements. 


3.1.  Emission  measurements 

All  measurements  were  conducted  in  a  non-intrusive  manner  using  optical  diagnostics  such  as 
emission  spectroscopy.  This  technique  is  not  spatially  resolved  but  since  the  plasma  presents  a  high 
degree  of  axisymmetry,  local  values  can  be  inferred  from  line-of-sight  measurements  along  several 
chords  of  the  plasma  circular  cross-section  through  the  Abel-invcrsion.  This  mathematical  transform 
imposes  several  constraints  on  the  optical  system  design.  First,  measurements  are  required  at  several 
lateral  locations,  preferably  every  millimeter  since  the  plasma  has  a  25  mm  radius.  Second,  the 
probed  volumes  must  be  constant  area  cylinders.  Hnaiiy,  to  allow  for  spatial  resolution,  the  probed 
volume  must  be  narrow,  ideally  pencil-like.  The  light  is  collected  by  a  50  cm  focal  length  UV  grade 
fused  silica  quartz  lens  placed  at  one  focal  length  from  the  jet  axis.  It  is  coupled  with  a  solid  angle 
limiting  iris  of  adjustable  circular  aperture  to  stop  the  beam  down  to  anywhere  above  f  25,  and  is  set 
on  a  stage  with  both  lateral  and  axial  translation  capability.  The  collimated  beam  is  then  refocussed 
by  a  15  cm  focal  length  UV  grade  fused  silica  quartz  lens  onto  the  monochromator  entrance  aperture 
which  acts  as  an  area  field  stop.  The  corresponding  system  magnification  is  0.3,  which  means  that 
the  imaged  plasma  area  is  3.3  times  the  monochromator  aperture  area.  With  a  typical  setting  of 
2  mm  for  the  gate  height  and  60  |im  for  the  slit  width,  a  spatial  resolution  of  0.2  mm  laterally  and 
6.6  mm  axially  is  achieved. 

The  light  then  goes  through  a  Stanford  Research  System  model  SR540  optical  chopper  into  a 
SPEX  model  1400-11  3/4  meter  scanning  monochromator.  The  choice  of  lenses  and  stop  fields 
ensures  proper  underfilling  of  the  7.5  entrance  f  number  monochromator.  The  1.200  lines/mm.  iOx 
10  cm  grating  is  blazed  at  3,000  A.  With  the  monochromator  used  in  single  pass,  the  linear  disper¬ 
sion  varies  between  11.1  .Vmm  at  2.000  A  and  10.4  A/mm  at  8.000  A.  Both  the  entrance  and  exit 
slits  are  continuously  adjustable  from  20  pm  to  three  millimeters.  The  dispersed  light  is  then  de¬ 
tected  by  a  Hamamatsu  model  R1104  head-on  photomultiplier  tube  with  a  high  photocathode  effi¬ 
ciency  spatial  uniformity  and  a  broad  quantum  efficiency  range  as  will  be  seen  in  Section  3.3.  Data 
are  acquired  using  a  Stanford  Research  Systems  model  SR510  lock-in  amplifier  and  then  transferred 
via  an  RS232  line  to  a  laboratory  computer  for  analysis. 

The  set-up  described  above  is  adequately  suited  for  emission  measurements  above  3.000  A. 
Measurements  in  the  ultra  violet  region  between  2,000  and  3,000  A  are  subject  to  specific  problems 
that  must  be  addressed  by  the  choice  of  different  optics.  Chromatic  aberration,  while  rather  miid  for 
most  lens  materials  in  the  visible  or  the  infra-red.  increases  dramatically  below  2,500  A  as  can  be 
seen  in  Fig.  3.2.  For  measurements  at  a  tlxed  waveiengm,  the  problem  can  oe  curumveniea  by 
adjusting  the  lens  position,  but  this  solution  is  not  practical  for  scans  over  large  wavelength  ranges. 
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Figure  3.1.  Experimental  set*up  for  emission  measurements. 

Among  the  possible  solutions  are  achromats.  but  they  are  difficult  to  build  for  wavelengths  be¬ 
low  2,500  A,  or  spherical  mirrors  in  lieu  of  lenses,  but  this  requires  major  changes  in  the  optical  sys¬ 
tem.  .A  simple  solution  is  to  minimize  the  chromatic  aberration  by  using  a  different  lens  matcriaL 
.Magnesium  fluoride  lenses  offer  a  dramadc  improvement  other  fused  silica  (Fig.  3.2)  while  stiil  al¬ 
lowing  for  a  broad  spectral  bandpass  (Fig.  3.3). 
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Fig.  3.2.  Chromatic  aberration  of  UV  grade  fused  silica  and  magnesium  fluoride.^ 

"he  iccoi  length  oi  a  lens  vanes  ’.viih  the  rciracuve  indc.\  'vnich  ;s  ;Lscil’  a  tuncuon  oi  ■'.aveiengih.  '’-.a 
vananons  are  very  well  approximated  by  the  thin  lens  formula;  f-./fx^  =  )/(n-v-l).  Refracuve  indices  tor  UV 

arade  fused  silica  and  magnesium  fluoride  are  taken  from  Oriel,  Optics  and  Filters,  p  12-5. 
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Fig.  3J.  Transmittivity  of  10  mm  thick  fused  siiica  and  3  mm  thick  magnesium  fluoride  windows.^ 

Another  problem  is  linked  to  the  use  of  the  UV  calibration  source.  The  deuterium  lamp,  dis¬ 
cussed  in  more  details  in  section  3.3,  has  a  very  small  circular  emitting  area  of  about  one  millimeter 
in  diameter.  To  insure  uniform  iiluminarion  of  the  entrance  slit,  the  magnification  of  the  system  has 
to  be  increased  in  order  to  at  least  double  the  image  size.  Since  this  increase  is  accompanied  by  a 
significant  signal  level  reduction,  a  compromise  must  be  reached.  For  these  reasons,  we  replaced  the 
collecting  lens  by  a  50  cm  magnesium  fluoride  lens,  and  the  imaging  one  by  a  100  cm  focal  length 
magnesium  fluoride  lens  for  all  UV  scans. 


3.2.  Radiation  source  strength  measurements 

Total  radiation  measurements  are  made  with  a  Scientech  model  36-0001  surface  absorbing  disk 
calorimeter.  The  detector  spectral  response  is  flat  between  2.500  and  25.000  A.  and  varies  by  about 
-1  %  between  2.000  and  2.500  A.  Self-calibration  is  provided  by  an  internal  heater  coil  to  an  accu¬ 
racy  of  about  3  %.  Lateral  measurements  are  performed  every  millimeter  using  a  translational  stage 
including  a  25  cm  focal  length  quanz  lens  used  in  2f/2f.  and  a  rectangular  area  stop  yielding  a  spatial 
resolution  of  1  mm  laterally  and  6.6  mm  axially.  The  6.6  mm  height  was  so  chosen  as  to  provide  the 
same  axial  spatial  averaging  as  for  temperature  measurements. 

Since  the  detector,  the  display  meter,  and  the  electrical  connexions  and  cables  are  very  sensitive 
to  the  RF  noise  generated  by  the  plasma  torch,  all  equipment  is  shielded  by  copper  Faraday  type 
cages.  However,  the  shielding  is  not  fully  efficient  in  removing  ail  the  noise.  .A  simple  solution  was 
found  with  the  introduction  of  a  low-pass  electrical  filter  between  the  analog  display  recorder  and  the 
digital  conversion  pon  of  the  lock-in  amplifier.  With  a  frequency  cut-off  at  40  Hz,  all  60  Hz  and 


-  Oriel.  Optics  and  Filters,  p  12-7.  8. 
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high  frequency  noise  are  eiiminaied  and  the  drift  rate  that  wouid  otherwise  be  sometimes  compara¬ 
ble.  over  several  minutes,  to  the  magnimde  of  the  measured  signal,  is  cut  down  to  less  than  a  percent 
of  the  measured  signal.  This  improvement  is  pardcuiarly  important  as  the  slow  3  s  response  time  of 
the  detector  causes  long  measurement  times. 


Otti  AequitU/an 
Computer 


Figure  i.4.  Experimental  set-up  for  radiation  source  strength  measurements. 

3.3.  Calibration 

Spectroscopic  temperature  measurements  require  the  knowledge  of  either  the  relative  or  the 
absolute  intensity  of  spectral  features.  The  measured  emission  spectra  need  also  to  be  corrected  for 
the  spectral  response  of  the  detection  system  which  includes  lenses,  minors,  monochromator  and 
photomuitiplier  tube.  The  detection  system  had  therefore  to  be  calibrated  in  intensity  between  2.000 
and  8.000  A.  To  cover  this  broad  spectral  range  with  maximum  accuracy,  the  combination  of  a  deu¬ 
terium  iamp  and  a  quartz  halogen  lamp  was  used.  Deuterium  lamps  cover  the  ultraviolet  range  2.000 
to  4.000  A  and  have  their  maximum  spectral  radiance  around  2.000  A  (see  Fig.  2.5).  Tungsten 
lamps  cover  the  rest  of  the  range  from  2.500  A  and  above  with  maximum  spectral  radiance  at  long 
■.vaveiengins  tsee  also  Fig.  3.5  ).  Tlie  overlapping  region  from  2.500  to  4.000  A  provides  a  check  tor 
the  quality  of  the  calibration. 
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Figure  3  Spectral  radiance  of  the  Optronic  tungsten  and  deuterium  standards. 

Both  lamps  used  in  this  work  are  standards  of  spectral  radiance  with  calibrations  traceable  to 
the  National  Insdmte  of  Standards  and  Technology  (NIST)  standards.  The  Optronic  model  UV- 
40IR  lamp  consists  of  a  40  Watt  deuterium  lamp  with  a  ultraviolet  transmitting  suprasil  window.  It 
is  calibrated  in  intervals  of  100  A  from  2,000  to  4,000  A  with  an  uncertainty  less  than  10%.  The  Op¬ 
tronic  model  5 50 A  lamp  is  a  modified  GE-I8A/T10/2P  tungsten  ribbon-filament  placed  behind  a 
sapphire  window  for  optimum  ultraviolet  transmission.  It  is  calibrated  from  2.500  to  8.000  A  in  in¬ 
tervals  of  100  A  and  has  a  quoted  uncertainty  of  less  than  3%. 

The  relative  spectral  response  of  the  detection  system  (opiicai  train  with  MgFj  lenses,  mono¬ 
chromator,  and  PMT)  is  presented  in  Fig.  3.6  where  the  plotted  curve  represents  the  normalized 
measured  output  signal  in  volts  per  unit  emitted  power.  Note  the  excellent  agreement  between  the 
two  lamps  m  the  overlapping  region  3,000-4.000  A. 

The  spectral  response  fall-off  at  the  2.000  A  end  is  due  to  the  photocathode  quantum  efficiency 
decrease,  to  the  decreasing  transmitdvity  of  the  two  lenses,  and,  mainly,  to  the  poor  reflectivity  of 
the  six  aluminum  silicon  monoxide  mirrors;  at  longer  wavelengths,  the  decrease  is  principally  caused 
by  the  poorer  quanmm  efficiency  of  the  photocathode  (Fig.  3.71.  The  two  peaks  around  5100  A  and 
6300  A  may  be  attributed  to  the  different  spectral  response  of  the  grating  to  the  parailely  and  per- 
penaicuiariy  poiarized  components  of  the  light. 
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Figure  3.6.  Relative  spectral  response  <  VoltAVatt)  of  the  detection  system  including  2  magnesium  lluoride 

lenses.  4  mirrors,  monochromator  and  photomultiplier  tube.  The  curve  is  normalized  bv  the  value 
at  5100  A. 

As  the  standards  of  spectral  radiance  are  calibrated  for  only  20  hours  of  use.  it  is  advantageous 
to  substitute  non-calibrated  tungsten  ribbon  lamps  such  as  the  GE-lSA/TlO/l  projection  lamps  for 
measurements  above  3,500  A.  In  this  case,  the  ribbon  temperature  is  measured  using  a  Pyro-Micro 
disappearing  filament,  single  color  pyrometer,  and  corrected  for  tungsten  emissivity  by  an  iterative 
procedure.  The  emissivity  of  tungsten  as  a  function  of  temperature  and  of  wavelength  was  measured 
by  several  investigators,  in  particuiar  Larrabee  (1959')  who  also  provided  curve  fits  between  1600 
and  2400  K  for  wavelengths  between  3,500  and  8,000  A.  When  performed  with  this  latter  method, 
the  relative  calibration  of  the  detection  system  above  3,500  A  is  found  to  be  in  excellent  agreement 
with  the  one  obtained  with  the  standard  of  spectral  radiance.  For  absolute  calibrations,  discrepancies 
do  not  exceed  10%. 
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Figure  3.7.  Relative  spectral  response  of  the  individual  components  of  the  detection  system. 

(a)  External  transmittivity  of  one  magnesium  fluoride  lens.^ 

(b)  Reflectance  of  one  AKSiO)  mirror.'* 

(c)  Radiant  sensitivity  of  the  R1104  photomultiplier  tube.^ 

(d)  All  4  mirrors  and  2  lenses  included  in  the  optical  train. 

(e)  Complete  detection  system:  4  mirrors,  2  lenses,  and  PMT. 


Oriel  Corporauon.  Optics  and  Filters,  p  i  2-8 
*  Oriel  Corporation,  Optics  and  Filters,  p  5-2 
'  Eamamaxsxi.PhotomuitiplierTuhes 
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4.  Thermodynamic  state  of  the  plasma. 


4.1.  Equilibrium 

Thermodynamic  characterization  of  a  piasma  involves  the  determinadon  of  the  number  density 
of  the  various  species  present  in  the  piasma  — atoms,  molecules,  ions,  and  electrons,  and  of  the 
distribution  of  particles  among  the  available  quantum  energy  states. 

In  the  limiting  case  of  a  plasma  in  Complete  Thermodynamic  Equilibrium  (CTE'),  ail 
microscopic  processes  affecting  the  distribution  of  the  species  are  in  detailed  balance.  In  that 
case,  all  plasma  particles  have  a  Maxwellian  velocity  distribution  characterized  by  a  common 
temperature  T,  and  radiative  emission  is  described  by  the  Planck  blackbody  function  at  the  same 
temperature  T.  The  thermodynamic  state  of  the  piasma  is  then  completely  determined  by  two 
parameters,  the  temperature  T  and  the  pressure.  With  the  notable  exception  of  the  solar  piasma 
below  the  photosphere,  such  a  situation  is  not  encountered  in  common  man-made  plasmas  the 
dimensions  of  which  being  too  small  to  prevent  radiation  escape. 

The  other  limiting  case  is  one  where  the  population  distribution  over  energy  states  cannot  be 
described  by  a  simple  disffibution  such  as  the  Boltzmann  distribution.  Instead,  each  level  needs  to 
be  determined  individually  which  is  practically  impossible  considering  that  the  number  of  available 
energy  states  is  inrinite. 

Many  common  applications  however  fall  in  between  these  two  extreme  cases  and  a  more 
complete  discussion  co'  found  in  Mitchner  and  Kruger  1 1973)  or  Bibeiman  et  al  (1987).  For 
atmospheric  pressure  plasmas  such  as  the  one  smdied  in  this  work,  the  frequency  of  collisions 
between  particles  is  high  enough  for  the  free  electrons  to  transfer  their  kinetic  energy  to  heavy 
ipeaes  emcientiy.  Such  a  piasma  is  referred  to  as  a  thermal  plasma,  and  the  kineuc  temperature 
of  the  free-eiectrons  is  close  to  the  kineuc  temperature  of  the  heavy  species,  at  least  in  the 
current-free  region  downstream  of  the  RF  coil  where  the  measurements  are  conducted. 
Moreover,  since  the  spacing  between  adjacent  rotational  levels  is  small  in  comparison  to  the 
average  kinetic  energy  associated  with  translational  motion,  and  since  the  energy  exchange 
between  translational  and  rotational  motion  occurs  very  rapidly,  molecules  tend  to  follow  an 
equilibrium  distribution  over  the  rotational  levels.  The  characteristic  temperature  of  the 
distribution,  termed  is  close  to  the  kinetic  (  translational)  temperature  T  of  the  gas. 

For  vibrational  levels,  the  simation  is  different  in  that  the  vibrational  energy  spacing  tends  to 
exceed  the  translational  energy  of  heavy  particles,  and  the  penod  of  molecular  vibrations  is  shon 
comoared  to  the  interaction  time  between  colliding  particles.  .\s  a  result,  vibrational-translational 
!  V-T)  energy  exchange  is  hindered  at  least  for  the  lowest  vibrational  levels  (for  higher  levels,  the 
period  of  molecular  vibrations  increases  and  the  V-T  exchange  increases  significantly).  Collisions 
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with  free  electrons,  on  the  otner  hand,  are  much  more  emcient.  and  eiectronic-vibranonai  energy 
exchange  tends  to  be  much  more  important  than  V-T  exchange.  In  addition,  vibrationai- 
vibranonai  energy  exchange  is  also  very  efficienL  As  a  result,  the  viorational  energy  levels  tend  to 
follow  a  quasi  equilibrium  distribution,  but  at  a  characteristic  temperature  T.j  close  to  the  free 
electron  temperature  T^.  Such  a  distribution  (Treanor  distribution)  holds  in  practice  for  the 
lowest  vibrational  levels  only  since  for  higher  lying  levels  the  rate  of  V-T  exchange  dominates, 
and  therefore  the  higher  levels  tend  to  be  closer  to  an  equil.ibrium  distribution  at  the  temperature 
of  the  heavy  praticies.  .As  will  be  seen  in  Section  5.  however,  only  a  few  vibrational  levels  aie 
significatly  populated  (typically  up  to  v'=4')  under  our  experimental  conditions  and  therfore  the 
vibrational  distribution  should  be  closer  to  the  free  electron  temperature. 

The  population  of  the  electronic  levels  is  governed  by  collisions  with  free  electrons,  and  by 
several  phenomena  such  as  radiation  escape  and  predissociation  that  cause  a  depletion  of  these 
ieveis.  Radiation  escape  arises  when  the  plasma  extent  is  finite;  m  that  case,  the  radiative  neid  is 
no  longer  blacbody  at  all  wavelengths.  Predissociation  is  a  molecular  process  by  which  the 
molecule  undergoes  spontaneous  (and  non-radiativei  disscxiation  because  of  a  crossing  with 
another  potential  curve.  .At  relatively  high  pressure,  and  atmospheric  pressure  in  particular,  the 
rates  of  coUisional  excitation  and  deexcitation  by  electron  impact  tend  to  be  much  higher  than  the 
depleting  radiative  and  predissociation  rates  (with  the  notable  exception  of  the  electronic  C  state 
of  NO  as  will  be  seen  in  Section  51.  As  a  result,  the  electronic  levels  tend  to  follow  a  Boltzmann 
distribution  at  the  temperature  of  the  free  electrons.  The  situation  is  however  complicated  in  such 
situations  as  ionizing  or  recombining  plasmas.  Since  high  lying  bound  electrons  undergo  little 
influence  from  the  nucleus,  they  tend  to  behave  in  the  same  way  as  free-eiectrons.  Therefore 
•hese  high  lying  levels  tend  to  oe  in  Saha  equilibrium  with  the  free  electrons  at  the  temperature  of 
the  free  electrons  but  not  inequiiibrium  with  the  ground  state.  Under  ionizing  conditions,  free 
electrons  are  produced  whereas  in  recombining  plasmas  they  undergo  depletion  through  ion 
recombination  processes.  Since  ionization  and  recombination  rates  are  finite,  a  defect  or  an 
excess  of  free  electrons  can  appear  because  of  inefficient  communication  between  the  low  and 
high  lying  bound  electronic  energy  levels.  In  the  case  of  a  recombining  plasma,  overpopulated 
high-lying  electronic  levels  can  be  caused  by  the  limited  rate  at  which  these  high-lying  levels 
cascade  down  to  lower  levels  (see  for  example  Owano  et  al.,  1990,  Gordon,  1992).  Situations 
such  as  those  depicted  in  Figure  4. 1  may  occur.  Only  in  the  particular  case  where  the  bound 
electrons  follow  a  Bolztmann  distribution  can  the  concept  of  the  electronic  temperature.  T^/,  be 
introduced. 

rnermai  equiiibrium  wiii  be  acnievea  u  ail  energy  leveis  can  be  uescnbed  by  a  common 
temperature  T.  If.  in  addition,  the  various  species  present  in  the  plasma  are  in  chemicai 
equilibrium,  the  plasma  is  completely  characterized  by  two  parameters,  the  temoerature  T  and  the 

:-I9 


pressure.  The  situation  is  anaiogous  to  the  case  or  CTE.  except  mat  the  concept  now  appiies  oniy 
locally,  and  only  to  the  species  distribution,  not  the  radiative  field.  This  case  is  commonly  referred 
to  as  Local  Thermodynamic  equilibrium  tLTEi. 


Figure  4.1.  Electronic  popuiation  distributions  in  equilibrium,  ionizing,  and  recombining  plasmas. 


4.2.  Thermal  equilibrium 

4.2.1.  Electronic  temperatures.  Population  distribution  of  the  electronic  levels. 

In  a  plasma  in  Local  Thermodynamic  Equilibrium  (LTEL  the  distribution  of  electronic  levels 
r'ollows  a  Boltzmann  statistic  characterized  by  a  single  temperature  T.  The  popuiation  of  any 
given  electronic  level  i.  with  energy  Ej  and  degeneracy  gj,  is  related  to  the  ground  level  popuiation 
:V)  ot  degeneracy  gQ  oy; 
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Taking  me  logariihm.  Eq.  4.1  becomes: 


,  f  n;  1  ,  f  no  ]  e. 

In  _.  =  in  —  j— — . 


Uoj 


(4.2) 


-A  plot  of  ]n{ni/gi)  versus  £,•  (Boltzmann  plot)  therefore  yields  a  straight  line  of  slope  -1/kT  from 
which  the  characteristic  temperature  of  the  distribution  can  be  readily  inferred. 

However,  as  stated  earlier,  this  ideal  situation  does  not  always  hold.  In  argon  plasmas  for 
instance,  the  upper  electronic  levels  may  be  overpopulated  with  respect  to  the  equilibrium 
distribution  when  the  electron  number  density  is  itself  elevated  (Gordon.  1992).  More  generally, 
depanures  from  LTE  occur  in  ionizing  or  recombining  plasmas.  In  our  case,  the  air  plasma 
ionizes  within  the  coil  region  and  recombines  as  it  cools  after  leaving  the  indu>,rion  zone.  Since  ail 
our  measurements  are  conducted  one  centimeter  downstream  of  the  nozzle  exit,  we  are  mainly 
concerned  with  recombination  effeas. 

When  the  plasma  cools,  the  free  electrons  tend  to  recombine  with  the  ions,  but  since 
recombination  rates  are  finite  the  process  is  not  instantaneous  and  an  excess  of  free  electrons  may 
appear.  The  outer  bound  electrons,  which  undergo  little  influence  from  the  nucleus  and  thus 
behave  similarly  to  the  free  electrons,  also  see  their  density  increase.  If  they  cannot  "cascade" 
down  fast  enough  into  the  lower  levels  of  the  particle,  the  high-  lying  electronic  levels  become 
overpopulated,  and  the  population  distribution  is  no  longer  Boltzmannian. 

In  order  to  investigate  a  possible  departure  from  equilibrium  in  the  electronic  levels,  two 
approaches  can  be  followed.  .A  first  method  (Sections  4.2.1.a  and  4.2.1. b)  consists  m  measuring 
;he  number  densities  of  several  electronic  levels  to  determine  whether  they  follow  a  Boltzmann 
distribution.  As  mentioned  earlier,  departure  from  a  Boltzmann  distribution  could  be  caused  by 
an  elevated  free-electron  number  density.  The  second  method  (Section  4.2. l.c)  therefore  consists 
in  measuring  the  free  electron  number  density  and  comparing  it  to  equilibrium  predictions. 


4.2.i.fl  LTE  Temperature 

The  LTE  temperature  is  defined  from  the  Boltzmann  relation  between  an  excited  electronic 
!evel  and  the  electronic  ground  state  of  the  particle: 
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Graphically,  it  is  obtained  on  a  Boltzmann  plot  from  the  slope  of  the  line  joining  the  ground  state 


to  the  eiecironic  ievei  under  consideranon.  For  a  piasma  in  LTE  this  temperature  is  ievei- 
independent  because  ail  levels  lie  on  a  straight  line.  Therefore,  u  the  piasma  is  in  LTE.  the  LTE 
temperature  does  represent  the  electronic  Ievei  distribudon.  On  the  other  hand,  if  the  electronic 
levels  do  not  follow  a  Boltzmann  distribution  with  respect  to  the  ground  state,  every  probed  level 
fields  a  different  LTE  temperamre.  In  this  case,  the  different  LTE  temperatures  have  no  physical 
meaning;  they  only  indicate  a  non-Boltzmann  distribution. 

The  number  density  of  an  excited  state  can  be  accessed  by  emission  spectroscopy  by  probing 
the  radiative  transitions  issued  from  that  state.  Transitions  from  atomic  electronic  levels  produce 
a  single  spectral  line  that  spans  a  narrow  spectral  range,  typically  a  few  angstroms.  .Atomic 
transitions  are  therefore  more  easily  accessed  than  molecular  transitions  that  span  very  large 
spectral  ranges,  typically  several  thousand  angstroms.  For  this  reason,  electronic  temperatures  are 
determined  in  the  present  work  from  atomic  transitions  only.  Not  all  atomic  transitions  can  be 
prooed  however,  they  have  to  lie  in  a  wavelength  range  adequate  for  the  detection  system  and 
they  must  be  easily  differentiable  against  the  underlying  background  due  to  continuum,  molecular 
bands,  or  second  order  light  These  restrictions  led  to  the  selection  of  the  five  atomic  transitions 
whose  characteristics  are  summarized  in  Table  4.1. 


Atom 

e/(cm-*) 

81 

£u  (cnr‘) 

8u 

X(A) 

Ay/  (10’  sec'*) 

flu 

Accuracy 

Oxygen 

73768 

5 

86629 

15 

7773.4 

3.40 

0.922 

10% 

Oxygen 

86629 

15 

102865 

25 

6157.3 

0.701 

0.0664 

10% 

Ninoeen 

83366 

6 

96752 

B 

7468.31 

1.61 

25% 

Nitroeen 

83319 

B 

96752 

B 

7442.30 

1.06 

0.088 

25% 

Nitroeen 

83286 

B 

96752 

4 

7423.64 

0.52 

0.086 

25% 

Table  4.1.  Electronic  transitions  considered  in  the  present  study  (Wiese  et  ai.  1966). 


The  total  intensity  I  of  a  spectral  line  is  related  to  the  number  density  of  a  transition  s  upper 
state  (emitting  level)  n.^  by: 


I  Wail  i_„  \ 

‘  — — 

[_cni"-srj  471 

where  the  Einstein  coefficient  and  the  upper  and  lower  energies  and  £;  are  the  constants 
indicated  in  Table  4. 1. 

The  procedure  for  determining  total  line  mtensities  is  as  follows.  First  the  monochromator  is 
lunea  to  me  waveiengtn  or  the  transition  unaer  consideration,  ia  our  case,  me  monochromator 
entrance  and  exit  slits  were  set  to  40  and  1 .000  pm  respectively,  which  correspond  to  special 
widths  of  0.44  and  1 1  X.  respectively.  The  large  exit  sht  width  is  necessary  to  minimize  line-wing 


(4.4) 
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cut-oif  and  to  capture  over  99%  or  the  line  energy  tsee  tor  instance  Owano.  i991a  piOl).  It  also 
enables  to  encompass  the  oxygen  hne  at  7773.4  A  which  is  in  fact  a  triplet  spread  over  3.44  A. 
The  total  intensity  emitted  by  the  plasma  is  then  measured  every  millimeter  along  plasma  chords, 
and  a  typical  lateral  intensity  profile  is  presented  in  Figure  4.2.  The  monochromator  is  then 
demned  by  a  few  angstroms  to  a  region  free  of  lines  where  the  measured  signal  corresponds  to 
background,  continuum,  or  second  order  light  Subtracting  this  signal  from  the  total  intensity 
produces  a  lateral  profile  of  line  intensity  only. 


Lateral  Position  (cm] 

Figure  4.2.  Lateral  emission  protiie  ot  the  oxygen  line  at  7.773  A 

This  profile,  in  unit  of  volts,  corresponds  to  the  line  intensity  averaged  along  plasma  chords. 
It  is  convened  to  radial  vanadons  using  a  mathemadcai  treatment  known  as  the  Abel-inversion. 
for  which  a  detailed  account  can  be  found  in  Appendix  1.  At  this  stage,  the  units  become  volt/m. 
Conversion  into  radiance  units  is  obtained  by  intensity  calibration  using  the  tungsten  lamp 
described  in  Chapter  3.  To  this  end.  the  calibrated  lamp  is  placed  at  the  focal  point  of  the 
detecting  lens  with  the  same  wavelength,  slit  settings,  and  PMT  voltage  as  for  line  intensity 
measurements.  The  procedure  \nelds  a  calibradon  factor,  in  Volt/( Watt/cm Vsri,  that  is  used  to 
transform  the  radial  intensity  profile  into  the  desired  radiance  units  of  Watt/cmVsr.  The  number 
density  of  the  emitdng  state  is  then  readily  obtaineo  from  £q.  4.4. 

The  remaining  parameter  needed  to  calculate  an  LTE  temperature  profile  from  Eq.  4.3  is  the 
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number  density  of  the  ground  state.  This  parameter  cannot  be  measured  by  emission 
spectroscopy — the  ground  state  does  not  emit  radiation,  nor  by  absorpdon  spectroscopy — the 
corresponding  electronic  transitions  are  in  vacuum  ultraviolet.  However,  a  measurement  is  not 
needed  since  the  ground  state  number  density  is  related  to  the  total  number  density  of  the 
species  s  (oxygen  or  nitrogen  here)  by  the  following  relation: 

_  nATiTE) 
go  QtiiTLTE) 

where  Q^i  stands  for  the  electronic  partition  function  of  the  atom  s.  Introducing  Xr,  the  mole 
fraction  of  species  s,  and  using  the  LTE  assumption. 

ns  =  XATlte)~ - ■  (4-6) 

IcTlte 

Subsntutmg  Eqs.  4.5  and  4.6  into  £q.  4.3; 


—  =  Xs(.TlTE)- 


I 


-e.xp  - 


(4.7) 


gu  lTE  QeiiT LTE)  "  LTE  J 

Eq.  4.7  is  an  implicit  equation  from  which  Tlte  ^^n  be  obtained  by  iteration.  Since  the  mole 
fraction  of  atomic  oxygen  or  nitrogen  varies  relatively  slowly  in  the  temperature  range  of  interest 
here  (see  Fig.  4.3),  the  behavior  of  Eq.  4.7  is  mainiy  controlled  by  the  exponential  termu  and 
therefore  convergence  is  reached  in  a  few  iterations,  typically  less  than  five. 


Temperaiure  [K] 


Figure  4.3.  Atomic  oxygen  mole  fraction  in  chemical  equilibrium  air. 


The  five  radial  LTE  temperature  profiles  thus  determined  are  presented  in  Fig.  4.4.  Each 
curve  has  been  truncated  at  the  radial  position  beyond  which  the  measured  PMT  signal  drops 
below  the  detection  limiL  The  very  good  agreement  between  the  five  temperature  profiles 
evidences  the  equilibrium  distribudon  of  the  electronic  levels.  However,  an  uncertainty  analysis 
and  an  estimate  of  the  effects  of  line  self-absorption  are  warranted  to  establish  the  validity  of  the 
temperature  profiles  as  accurate  descriptions  of  electronic  level  distribution. 


Figure  4.4.  LTE  temperature  profiles. 


Uncertainty  analysis 


Special  attention  was  paid  to  determine  the  error  bars  on  the  temperature  profiles.  For 
that  purpose,  a  computer  program  was  written  to  simultaneously  perform  the  data  reduction  and 
an  uncertainty  analysis  that  accounts  for  the  standard  deviation  on  the  measurements,  the 
uncertainty  propagated  by  the  fitting  and  the  Abel-inversion  procedures,  the  uncertainty  on  the 
absolute  calibration  in^pically  5%),  and  the  uncertainty  on  the  Einstein  coefficients  (see  Table 
-i.l).  The  error  bars  thus  determined  are  presented  in  Figure  4.5  for  the  five  atomic  lines. 
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These  error  bars,  which  are  typically  on  the  order  of  100  K  or  1.5%,  may  appear  at  first 
sight  surprisingly  small.  The  following  simple  analysis,  however,  should  shed  some  light  on  the 
reasons  for  such  small  uncertainties. 

Taking  the  logarithm  of  Eq.  4.7  and  differentiating,  we  obtain: 
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The  electronic  partition  functions  of  atomic  oxygen  and  nitrogen,  respectively  equal  to 

(22  =  5  +  +  o(^-i3.ooorr)^  (4  9^ 


(22  =  4+o(e-^8.ooo/r)^  ,4  jO) 

are  practically  constant  over  the  temperature  range  of  this  study  and  can  be  considered  equal  to  9 
and  4.  respectively,  with  better  than  1%  accuracy  up  to  10.000  K.  .\lso.  as  may  be  seen  from 
Fig.  4.3.  the  oxygen  and  nitrogen  mole  fractions  vary  slowly  with  temperature.  Hence, 
neglecting  \^X^/  X and  | A  Q,,  /  QJ : 
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Considering  for  instance  the  oxygen  line  at  7,773  A,  the  uncertainty  at  an  average  temperature  of 
7.000  K  is  therefore  given  by: 


^Tlte\_  1 

^{njgj\ 

Tlte  17.8-1-1 

njg^  1 

=  0.051 
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Consequently,  at  7,000  K,  it  is  sufficient  to  know  within  30%  accuracy  to  obtain  no  more  than 
1.5%  (or  lOOK)  uncertainty  on  the  LTE  temperamre.  At  this  rate,  a  difference  of  several 
hundred  degrees  between  the  different  LTE  temperatures  would  mean  a  non-Boltzmannian 
distribution.  This  is  not  the  case  however  for  any  of  the  five  temperature  profiles  since  they  are 
typically  within  a  few  ten  degrees  of  each  other. 


Line  self-absorption 


As  it  crosses  its  way  out  of  the  plasma,  the  light  emitted  inside  the  plasma  undergoes  some 
absorption.  Since  the  absorption  coefficient  is  proportional  to  the  line  oscillator  strength  or 
equivalently  to  the  Einstein  coefficient.  ‘  the  values  listed  in  Table  4. 1  suggest  that  the  oxygen  line 
at  7,773  A  is  subject  to  the  most  intense  seif-absorption  of  the  5  lines.  For  a  plasma  in  LTE,  the 
spectral  absorption  coefficient  has  the  following  expression: 
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where  is  the  number  density  of  the  lower  level  of  the  transition,  the  electronic  oscillator 
strength,  and  <I>(A,)  the  line  shape  factor  that  includes  Doppler,  natural.  Stark,  van  der  Waais.  and 
resonance  broadening  mechanisms.  While  Doppler  broadening  results  in  a  Gaussian  lineshape: 
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'  The  Einstein  coefficient  A..;  [scc'M  is  related  to  the  dimensionless  absoronon  oscillator  streneth  by: 
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the  four  omer  mechanisms  produce  a  Lorentzian  pronie: 


(4.15) 


We  have  neelected  in  the  above  equation  the  frequency  shift  that  usually  accompanies  Lorentzian 
mechanisms  for  it  is  much  smallf»r  than  the  linewidth  of  the  transitions  considered  in  our  study.  It 
should  be  noted  also  that  the  total  Lorentzian  half  width  is  the  sum  of  ail  Lorentzian  widths.  i.e. 
AA.L  =  ^.uuurai+^Stark-^^Resommce^^:ander  Waais-  The  pronie  due  to  both  broadening  types 
is  called  the  Voigt  profile: 
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;sVln2  ,  and £  =  vln2^  -^. 
AXd  '  AXd 
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It  is  therefore  important  to  determine  which  broadening  mechanisms  have  an  influence  on  the 
total  line  broadening.  This  and  expressions  for  the  half  widths  at  half  maximum  are  presented 
below.  All  terms  in  the  equations  are  in  SI  units,  unless  explicitly  specified  otherwise. 

Doppler  broadening:  If  the  photon  emitted  during  the  transition  has  a  frequency  Vg  in  the 
reference  frame  of  the  emitting  molecule,  an  observer  placed  in  the  laboratory  frame  will  see  the 
molecules  coming  in  its  direction  with  a  velocity  v  emit  a  photon  of  frequency  v  =  vo(/+v/c).  If 
the  molecules  have  a  Boltzmann  veiocitv  distribution  characterized  by  the  heavy  particle 
:emperature  the  resulting  Doppler  line  is  Gaussian  with  a  half  width  at  naif  maximum  given  by: 


Natural  broadening  is  due  to  the  Heisenberg  uncertainty  principle  applied  to  the  enei^es  of 
the  initial  and  the  final  state  of  the  transition.  This  results  in  an  uncertainty  on  the  transition 
frequency  or  wavelength  which  translates  into  a  Lorentzian  line  profile  of  half  width  at  half 
maximum: 


^  =  3.58x10"^ 
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where  the  two  summation  terms  are  the  inverses  of  the  transition  s  upper  and  lower  level  lifetimes. 
Natural  broadening,  which  is  generally  negligible  in  comparison  with  the  other  mechanisms. 


represents  a  lower  bound  to  the  spectral  linewidth.  In  the  case  of  the  oxygen  line  at  7.773  A  the 
natural  half  width  at  half  maximum  was  calculated  using  the  constants  of  Wiese  et  al  (1966)  and 
found  to  be  ^saurai  ~  3.4x10'^  A. 

The  next  three  broadening  mechanisms,  resonance,  van  der  Waais.  and  Stark,  are  respectively 
due  to  collisions  with  "like"  neutrals,  "unlike"  neutrals,  and  charged  particles.  A  brief  descripdon 
follows  along  with  numerical  expressions  and  applications  to  the  oxygen  line. 

Resonance  broadening  is  caused  by  collisions  between  "like"  particles  (e.g.  two  oxygen 
atoms  )  where  the  perturber's  initiai  state  is  connected  by  an  allowed  transition  to  either  the  upper 
or  lower  state  of  the  radiative  transition  under  consideration.  Typically,  the  three  perturbing 
transitions  that  need  to  be  considered  are  g-^l.  g-^u.  and  /-»w.  where  g  stands  for  the  ground 
state  of  the  particles,  and  I  and  u  for  the  lower  and  upper  states  of  the  radiative  transition,  as 
usual.  Using  the  expression  given  by  Griem  ( 1964,  p  97),  we  obtain: 
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For  the  oxygen  line  at  7,773  A,  the  constants  listed  by  Wiese  et  al  (1966)  for  the  relevant 
transitions  lead  to: 


[  A  ]  =  (2.76  X  lO'^  +  6.851exp(-106.226  /  T,  ))Xo7"[atmjr^,  (4.21) 

which  has  a  maximum  of  3.5x10'^  .4  at  10,000  K.  In  other  words,  resonance  broadening  is 
uompieteiy  negligible  for  this  oxygen  transition — indeed  only  the  lower  state  is  connected  to  the 
ground  state  by  an  allowed  transition,  and  only  very  weakly  so  since =  5.4x10'^. 


Van  der  Wools  broadening  is  caused  by  collisions  with  neutral  pemirbers  that  do  not  share  a 
resonant  transition  with  the  radiating  particle.  Griem  (1964,  p  99)  gives  the  following  expression 
for  a  radiator  r  colliding  with  a  pemirber  p: 
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In  the  above  expressions.  Ep  is  the  energy  of  the  first  excited  state  of  the  perturber  connected  with 
its  ground  state  by  an  allowed  transition.  £//  and  are  the  ionization  potential  of  the  hydrogen 
atom  and  of  the  radiating  atom,  respectively,  £q  is  the  excitation  potential  of  the  upper  state  of 
the  line.  is  its  orbital  quanram  number,  z  is  the  number  of  effective  chaises,  Np  is  the  number 
density  of  the  perturbing  species,  and  v  is  the  relative  speed  of  the  radiating  atom  and  the 
perturber.  The  relative  velocity  term  can  be  expressed  as  a  function  of  the  mean  speed: 


vl'/  =  {4/7c)''‘°  r(9/5)  (v, j'''  =  0.98  (v,,f '  =  0.98  (^8£r/7tm; (4.24) 


where  ml.  is  the  reduced  mass  of  the  radiator  and  the  perturber.  Summing  over  ail  perturbers 
present  in  the  plasma  and  introducing  the  mole  fraction  Xp  of  perturber  p,  Eq.  4.22  becomes: 
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The  summation  in  this  expression  runs  over  aU  neutral  species  present  in  the  plasma.  For  an  air 
plasma,  the  five  perturbers  that  must  be  considered  are  O,  N,  N2,  O2,  and  NO,  which  together 
constimte  98%  of  the  chemicai  equilibrium  composition  for  temperatures  up  to  10, (XX)  K.  We 
computed  the  mole  fractions  of  these  five  components  up  to  lO.CXX)  K  and  combined  these  values 
with  the  EpS  and  that  are  indicated  in  Table  4.2.  Note  that  the  reduced  masses  indicated  in 


the  table  are  for  the  case  where  the  radiator  is  an  oxygen  atom.  When  a  different  radiator  is 
considered,  a  new  set  must  be  considered. 


Perturber 


[g/mole] 


O 


8.00 


O'. 


7.47 


10.67 


N'. 


10.18 


Transition  issued  from  the  Hrst  excited  state 
opdcailv  connected  to  the  ground  state. 


£;,[eVl 


3s0^p0 


9.5 


0.088 


4p  _^4s0 


A3Z/ 


X^Z  -  (Herzberg  1) 


a'^Zy'  ^  X'Zjj*  (WilUdnson-Mulliken) 


10.3 


0.085 


4.4 


0.150 


8.5 


0.090 


NO 


10.43 


A^Z"''  ^  X^n  (Gamma  ) 


5.5 


0.127 


Table  4.2.  Constants  needed  in  Eq.  4.25  when  the  radiator  is  an  oxygen  atom.  For  other  radiators,  the 
reduced  mass  must  of  course  be  modified  as  warranted. 
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The  summation  lenn  in  Eq.  4.25  was  calculated  up  to  10.000  K  for  several  radiators,  namely  O. 
N.  Oo.  Ni.  and  NO. 


As  can  be  seen  in  Fig.  4.6.  this  term  is  almost  constant  over  the  entire  temperature  range  for 
either  O,  N,  Ni  Oo  or  NO  as  a  radiator.  An  approximate  linear  fit  of  the  curves  shown  in  Fig.  4.6 
was  penbrmed  for  each  of  the  radiating  species  previously  mentioned.  For  au"  in  chemical 
equilibrium  at  temperatures  up  to  10.000  K.  it  is  therefore  possible  to  give  anaiyticai  formulas  for 
the  van  der  Waals  broadening  of  O,  N,  Oo,  N2,  and  NO  lines: 

^vanaerWaaU  =  4.0  XlO'^ )’ 7^ f  C, (T)  (4.26) 

where  depends  only  on  the  radiating  species  and  on  the  transition  under  consideration,  and 
where  CJT),  defined  as 


CAT)  = 


1 
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(4.27) 


is  a  function  that  only  depends  on  temperature  (although  very  weakly)  and  on  the  radiating 
species.  The  C/T)  functions  for  O.  N,  N2,  O2,  and  NO  are  summarized  in  Table  4.3. 
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Radiator  r  !  C/T)  [(eVV*^  *(g/moici'°-^1 


OorN 

0.10-1.9x10-^ 

N2,  Ot  or  NO 

0.08-1.3x10-6T 

Table  4.3.  Van  der  Waais  broadeniiig  coefficteDts  for  O.  N,  N2.  02>  and  NO  lines  in  a  chemical  equilibrinm 
air  plasma.  These  fits  are  within  better  than  10%  of  the  actual  values  for  all  temperatnres  up  to 
10,000  K. 


We  can  now  detennine  an  expression  for  the  van  der  Waais  broadening  of  the  oxygen  atom 
transition  at  7,773  A.  Recalling  that  for  oxygen  atoms  the  effective  charge  number  2  is  equal  to  1 
and  that  the  ionization  potential  is  13.614  eV,  that  the  excitation  potential  of  the  upper  state  of 
the  line  at  7,773  A  is  10.75  eV  and  its  orbital  quanmm  number  is  1.  and  finally  that  Efj  = 

1 3.6  eV.  the  matrix  element  of  Eq.  4.24  can  be  calculated: 

We  finally  obtain  for  the  van  der  Waais  broadening  of  the  oxygen  line  at  7,773  A; 

^.an^^waais  =  (l  1.2 -2.12x10-*  7)r^V[mm]  (4.29) 

Van  der  Waais  broadening  therefore  decreases  with  temperamre,  which  is  not  surprising  since 
although  the  collision  frequency  increases  with  temperature,  the  density  decreases  faster. 

Stark  broadening  is  due  to  the  Coulomb  interaction  between  the  outer  electron  of  the 
radiator  and  the  charged  particles  present  in  the  plasma.  Both  ions  and  electrons  induce 
broadening  but  because  of  their  comparatively  higher  rciative  velocities  electrons  are  responsible 
for  most  of  it.  For  this  reason,  we  consider  only  electron  Stark  broadening  which  is  a  function  of 
the  free  electron  number  density  and  a  weak  function  of  the  electron  temperature.  Griem  ( 1964'! 
gives  the  following  expression  for  electron  Stark  broadening: 

(4.30) 

where  wiTg)  is  a  tabulated  function  in  units  of  A  that  depends  on  the  radiator  and  on  the 
transition.  Griem  provides  numerical  values  for  the  most  important  transitions  of  a  number  of 
atoms  (Griem,  1964.  Table  4.5,  pp.  454  j5^.).  For  the  oxygen  line  at  7.773  A,  wiTg)  varies 
between  1.99x10-2  A  at  2.500  K.  2.48x10-2  A  at  5.000  K.  and  3.27x10-2  A  at  10,000  K  (ib.. 
p  4681. 
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(4.28) 
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Broadening  widths  for  the  mechanisms  presented  above  are  plotted  as  a  funcdon  or 
temperature  up  to  10,000  K  in  Fig.  4.7  where  it  appears  that  Doppler  and  Van  der  Waais 
broadening  are  the  dominant  mechanisms  for  the  oxygen  line.  A  Voigt  profile  must  therefore  be 
computed  and  subsdmted  in  Equadon  4. 13  to  determine  the  spectral  absorption  coefficient  as  a 
function  of  temperature.  Since  no  analytical  expression  exists  for  the  Voigt  profile,  the  numerical 
program  proposed  by  Drayson  ( 1976)  was  used  to  compute  Voigt  line  profiles. 


We  can  now  return  to  the  problem  of  interest  which  is  to  correct  the  measured  temperature  for 
line  self  absorpnon.  The  procedure  is  as  follows:  We  first  use  the  measured  lateral  intensity 
profile  and  assume  that  there  is  no  absorption  to  determine  an  "opdcally  thin"  temperature  profile 
tfom  the  line  under  consideration.  We  then  divide  the  plasma  in  25  equally  spaced  (1  mm) 
concentric  rings  in  which  thermodynamic  properties  are  nearly  constant  and  defined  by  the 
optically  thin  temperature  profile. 

With  the  conventions  of  Fig.  4.8.  the  radiative  transport  equation  for  a  given  ring  «  at  a  given 
lateral  location  z  reads: 

/,(X)  =  S,(rj[l-exp(-i<:,(X)/,(z))j  +  /,(X)exp(-i«:,(X)/,(2))  (4.31) 

where  k.^  is  the  absorption  coefficient  corrected  for  stimulated  emission,  which  is  related  to  the 
absorption  coefficient  by  the  following  relation  for  a  plasma  in  LTE: 
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(4.32) 


and  where  ByjT)  is  the  Planck  blackbody  function: 


Sx(7;)  = 


_ C _ 


with  Cl  =  1.19x10^  Wptrrcm'^sr^  and  =  1.438x10^  pmK. 


(4.33) 


Flgare  4.8.  Defining  elements  for  the  radiative  transport  equation  Eq.  4  Jl. 

In  Eq.  4.31,  the  first  term  corresponds  to  the  net  light  emitted  by  the  plasma  in  the  n^**  layer, 
and  the  second  term  corresponds  to  the  amount  of  light  incident  on  the  n^  layer  that  got 
transmitted  by  that  layer.  When  ail  layers  are  included,  a  spectral  profile  is  obtained  for  the  given 
lateral  position  z.  .-Mong  a  plasma  diameter  (r=0),  one  obtains  the  spectral  profile  shown  in 
Fig.  4.9.  At  the  exit  of  the  last  layer,  the  intensity  profile  is  spectrally  integrated  to  detenxdne  the 
total  line  intensity  corrected  for  seif  absorption. 

We  can  thereby  reconstruct  a  lateral  profile  of  net  emitted  intensity,  and  compare  it  to  what  it 
would  be  if  the  plasma  were  optically  thin.  This  procedure  was  numerically  performed  and  T|j(z), 
the  resulting  ratio  between  the  profile  including  spectral  absorption  and  the  optically  thin  one.  is 
plotted  in  Hg.  4. 10  for  the  oxygen  triplet  at  7,773  A.  We  then  multiply  the  measured  lateral 
intensity  profile  by  Tj^fz)  and  recompute  a  new  temperature  profile  using  the  same  method  as 
previously  described.  This  procedure  is  then  iterated,  that  is  a  new  profile  of  seif  absorbed  over 
optically  thin  intensity  Ti2(z)  is  calculated  to  more  accurately  correct  the  original  measured  lateral 
intensity  profile,  and  again  a  new  temperature  profile  is  calculated,  .\fter  a  few  iterations,  the 
temperature  profiles  do  not  change  by  more  than  a  fraction  of  a  degree  Kelvin.  At  this  point,  the 


2-34 


resulting  temperature  pronie  is  corrected  for  seif  absorpnon.  it  should  be  noted  here  that 
convergence  to  less  than  a  few  degrees  is  typically  obtained  with  only  one  iteration. 
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Fignre  4.9.  Numerical  prediction  of  the  spectral  emission  of  the  oxygen  triplet  at  7,773  A  along  a  diameter  of 
the  plasma  (z=0). 
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Figure  4.10.  Ratio  between  the  optically  thin  and  the  self-absorbed  integrated  emission  of  the  oxygen  triplet 
as  a  function  of  lateral  position  z. 


The  LTE  temperature  profile  of  the  oxygen  line  at  7,773  A  corrected  for  seif  absorption  is 
plotted  in  Fia.  4. 1 1  along  with  the  optically-thin  profile.  .As  can  be  seen,  the  difference  is  at  most 
100  K.  which  is  no  more  than  the  uncertainty  due  to  Einstein  coefficients  and  experimental  errors. 
A  qmilar  treatment  was  also  applied  to  the  four  other  temperature  profiles  and  the  corrections 
never  exceed  80  K. 

The  conclusions  reached  earlier  regarding  the  LTE  distribution  of  the  electronic  levels  are 
therefore  unchanged  after  an  uncertainty  analysis  and  corrections  for  line  self-absorpnon. 


Rgure  4.11.  Self  absorption  correction  of  the  LTE  temperature  from  the  7,773  A  oxygen  tripleL 


4.2.1. b  Boltzmann  temperature 

Another  way  to  examine  possible  departures  from  LTE  is  developed  in  this  section  with  the 
inttoduedon  of  the  so-called  Boltzmann  temperamre.  The  Boltzmann  temperamre  is  defined 
by  the  distribution  of  excited  electronic  levels.  It  is  obtained  from  the  Boltzmann  relation  between 
the  number  densities  of  at  least  two  excited  levels  of  a  same  atom. 

M  =  £i£expf-  .  (4.34) 

.Sul  Su2  ^  ^Tb  . 

Two  transitions  originating  from  different  upper  states  ul  and  u2  are  therefore  needed.  While 
the  oxygen  tines  of  the  previous  secuon  may  be  used,  the  nitrogen  unes  ao  not  qualify  because 
they  have  a  common  upper  state.  Although  LTE  and  Boltzmann  temperatures  are  not 
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independent  (Eq.  4.34  is  the  rano  of  Eq.  4.3  wntten  for  two  electronic  levels),  the  Boltzmann 
temperature  offers  several  advantages.  First,  it  has  a  physical  meaning,  because  even  when  the 
plasma  is  not  in  LTE,  Tq  represents  the  kinenc  temperature  of  the  free  electrons  if  the  probed 
levels  are  close  to  the  ionization  limiL  This  stems  from  the  fact  that  high-lying  electronic  levels 
tend  to  be  in  equilibrium  with  the  free  electrons.  Second,  in  contrast  to  the  LTE  temperature,  no 
assumption  needs  to  be  made  about  the  ground  state  number  density  since  it  is  no  longer  involved. 
Finally,  the  Boltzmann  temperature  is  more  sensitive  to  a  possible  non-equilibrium  distribution  of 
the  electronic  levels.  However,  the  latter  point  also  implies  that  the  Boltzmann  temperature 
cannot  be  determined  as  accurately  as  the  LTE  temperature  because  it  involves  the  Finwpin 
coefficients  of  two  transitions,  and  these  are  not  yet  known  very  precisely  (10%  for  the  two 


0  0.5  1,0  1.5  2.0  2.5 

Radial  position  (cm) 

Figure  4.12.  Boltzmann  temperature  determined  from  the  relative  intensity  of  the  two  oxygen  lines. 


The  following  uncenainty  analysis  gives  an  idea  about  the  precision  to  expect  for  the 
Boltzmann  temperature  based  on  two  lines.  Defining  Eq.  4.34  may  be  rearranged  as: 

-!-  =  — - — Inf-^l  (4.35^ 

Tb  £uz-Zu\  yriui ) 

Taking  the  logarithm  and  differentiating,  this  equation  becomes; 

'-r,  'aU:,;,,-.:)! 

£«2  £ut  j  flux! flu2  j 

The  energy  difference  between  the  two  oxygen  levels  being  2.03  eV.  the  uncenainty  at 
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1 4.36 1 
6.000  K  is 
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then 


ATfli, 

tb  r 


0.251 


nil/ nil  I 


(4.37) 


Since  the  number  density  ratio  cannot  be  known  with  better  than  20%  accuracy  due  to 
uncertainties  on  Einstein  coefficients  alone,  the  Boltzmann  temperature  cannot  be  determined  to 
an  accuracy  better  than  5%.  Although  it  can  be  argued  that  the  ratio  of  Einstein  coefficients 
should  be  more  accurate  than  their  absolute  magnitude,  the  error  bars  indicated  in  Fig.  4.12  do 
include  the  effect  of  10%  uncertainty  on  these  coefficients  in  addition  to  all  the  other  uncertainty 
causes  that  were  mentioned  previously  in  the  LTE  temperature  determination  procedure.  Also 
plotted  in  Fig.  4.12  are  the  LTE  temperatures  from  the  oxygen  lines  at  7.773  A  and  at  6,157  A. 
The  comparison  shows  that  the  Boltzmann  temperature  is  within  e.xpcrimentai  uncertainty  from 
the  LTE  temperatures  discussed  in  the  previous  section. 

The  last  question  to  be  addressed  here  concerns  the  effect  of  seif  absorption  on  the  Boltzmann 
temperature  profile.  The  two  oxygen  line  intensity  profiles  corrected  for  seif  absorption  lead  to 
the  Boltzmann  temperature  presented  in  Fig.  4. 13  along  with  the  corrected  LTE  temperature  firom 
the  oxygen  lines  at  7,773  A  and  at  6,157  A.  This  corrected  Boltzmaim  profile  is  also  within 
experimental  uncertainty  of  the  LTE  temperatures. 

In  conclusion,  the  Boltzmaim  temperature  also  suppon  a  Boltzmann  distribution  of  the 
electronic  levels. 


Figure  4.13.  Boitzmann  temperature  determined  from  the  relative  intensity  of  the  two  oxygen  lines 
corrected  for  self  absorption. 
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J.2.I.C.  Electron  number  density 


A  more  direct  way  to  determine  whether  the  outer  electronic  levels  are  overpopulated  due  to 
an  elevated  free  electron  number  density  is  precisely  to  measure  the  free  electron  number  density. 
To  that  purpose,  a  small  amount  (1.9  1/min)  of  hydrogen  Hi  was  mixed  with  the  95  1/min  air  flow 
and  injected  into  the  torch.  .At  the  operating  temperature  conditions,  H2  dissociates  almost 
completely  and  several  atomic  hydrogen  lines  from  the  Balmer  series  appear  in  the  visible,  in 
particular  the  transition  at  4861.32  A.  The  Hp  line  is  Stark-broadened  due  to  collisions 
between  the  hydrogen  atoms  and  the  charged  particles  which  are  present  in  this  temperature 
range. 

Stark  broadening  of  the  Hp  transition  has  been  widely  smdied  (Griem,  1964,  Vidal  et  al, 
1973).  It  is  a  function  of  the  free  electron  concentration  and  a  weak  function  of  temperature.  A 
fit  of  the  values  given  by  Vidal  et  al  (1973)  gives  for  the  Stark  Half  Width  at  Half  Maximum: 

,  f  0.647 

A]=1.740xl0-‘°(n,[cm--‘J)  (4.38) 

This  fit  is  valid  for  electron  densities  between  10*^  and  10*^  cm*^  but  it  does  not  take  into 
account  the  temperature  dependence  since  the  correction  would  not  exceed  a  factor  20%  over  the 
temperature  range  5,0(X)  to  40,0(X)  K. 

The  other  broadening  mechanisms  are  also  examined  following  the  same  method  as  in  Section 
4.2.1.b.  Half  Widths  at  Half  Maximum  for  the  Hj3  line  in  an  air  plasma  are  summarized  in  Table 
4.4  and  plotted  as  a  function  of  the  electron  number  density  in  Figure  4. 14.  The  conversion  from 
temperature  to  electron  number  density  was  done  assuming  a  chemical  equilibrium  concentration 
of  electrons.  Over  the  temperature  range  of  interest  here  (5,000  to  7.500  K),  resonance 
broaaemng  is  negligible  since  it  does  not  exceed  0.005  A  for  a  hydrogen  mole  fracnon  of  4%. 


^resonance 

‘^^an  der  Wools 

^Doppler 

(24.5-6.00x  10“*  T)^^5^ 

1.74x10"^  Vf 

Table  4.4.  Half  Widths  at  Half  Maximom  (in  A)  for  the  hydrogen  beta  line  at  4861 J2  A 


The  monochromator  entrance  and  exit  slits  were  both  set  to  40  qm  which  should  correspond 
to  a  Gaussian  slit  function  of  half  width  at  half  maximum  equal  to  0.22  A.  These  very  small 
Aidths  are  necessary  to  introduce  as  little  instrumental  broadening  as  possible.  However,  as  can 
be  seen  in  Fig.  4,14,  the  resulting  theoretical  instrumental  broadening  is  already  larger  than  the 
Doppler  broadening  which  is  the  only  other  mechanism  that  gives  a  Gaussian  component  to  the 
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Voigt  proriie.  It  is  therefore  cnticai  to  determine  the  insmimeniai  broaaemng  quite  accurately. 
The  quoted  instrumental  broadening  is  not  very  accurate  because  the  physical  slit  settings  are 
close  to  the  narrowest  achievable  width  for  the  monochromator.  For  that  reason,  the  actual  slit 
function  was  measured  with  a  narrow  spectral  source  (HeNe  at  6,328A);  it  was  found  to  be  nearly 
Gaussian  with  a  half  width  at  half  maximum  of  0.3  ±0.05  A.  An  insmimeniai  broadening  of  0.3  A 
was  therefore  taken  into  account  to  determine  the  relation  between  the  half  width  at  half 
maximum  and  the  electron  number  density.  The  resulting  Voigt  HWHM  is  shown  in  Fig.  4. 14. 
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Figure  4.14.  Broadening  widths  for  the  hydrogen  beta  line  at  4861.32  .4. 


The  Hfj  line-shape  was  measured  every  mm  over  the  radius  of  the  plasma,  yielding  25  pronles 
such  as  the  one  shown  in  Fig.  4.15.  The  distortion  is  due  to  the  underlying  continuum  which 
essentially  comprises  bands  of  the  second  positive  system  of  N2.  The  spectrum  was  therefore  also 
measured  in  pure  air  to  determine  the  baseline.  When  the  baseline  is  subtracted,  a  symmetric  line- 
shape  is  recovered.  The  line  profiles  were  then  Abel-invened.  and  the  electron  number  density 
determined  from  their  half  width  at  half  maximum. 

Figure  4.16  shows  the  measured  electron  number  density  profile  out  to  a  radius  of  about 
1.5  cm.  Beyond  that  point,  the  Hp  line  intensity  decreases  below  the  detection  limit  Also  shown 
:n  Figure  4.16  is  the  chemical  equilibrium  electron  de.nsity  which  was  calculated  from  the  LTE 
temperature  of  the  oxygen  line  at  7.773  A  corrected  for  self-absorption.  The  close  agreement 
between  the  two  profiles  shows  that  the  electron  number  density  is  not  elevated  with  respect  to  its 
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equilibrium  value.  This  election  number  density  measurement  therefore  adds  further  suppon  to  a 
Boltzmann  distribution  of  the  electronic  levels. 
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Figure  4.15.  Hp  line  profile  reconstruction  obtained  by  subtracting  the  baseline  signal  due  mainly  to  N2 
from  the  total  measured  intensity. 
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422.  Rotational  temperature 


4.22.0.  N2*(0,0)  band 

Rotational  temperatures  can  be  measured  by  resolving  rotational  lines  and  inferring  a 
temperature  from  molecular  Boltzmann  plots  (Tourin,  1966).  We  have  applied  this  technique  to 
the  rotational  lines  of  the  (0,0)  vibrational  band  of  the  Nj  first  negative  system.  However,  even 
with  many  rotational  lines  (23  in  our  case),  the  temperature  profile  obtained  through  Abel- 
inversion  is  extremely  sensitive  to  the  uncertainty  in  the  measurements  because  the  energy  spacing 
between  rotational  levels  is  very  small.  Consequently,  we  only  determined  a  temperature  firom  a 
rotational  spectrum  measured  along  the  diameter  of  the  plasma.  This  temperature  corresponds 
nevertheless  to  an  average  rotational  temperature  over  the  central  part  of  the  plasma  where  N2 
has  its  maximum  emission.  As  shown  in  Figure  4.17,  this  radially  averaged  temperature  lies 
within  the  LTE  and  Boltzmann  temperatures,  and  is  therefore  consistent  with  rotational  equilib¬ 
rium. 
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4.2J2.b.  NO  y(OJ)  band 

A  rotational  temperature  was  also  obtained  from  the  Abel-inverted  shape  of  the  (0,1)  band  of 
NO  Gamma  following  the  work  by  Gomes  et  al  (1992).  This  technique  is  currently  limited  to  re¬ 
gions  below  6,500  K,  since  above  that  temperature  NO  Delta  and  Epsilon  bands  interfere  with 
NO  Gamma.  The  technique  could  be  extended  above  6,500  K  if  the  absolute  intensity  of  these 
superimposed  bands  were  known  with  good  accuracy,  but  this  is  not  yet  the  case.  Again,  good 
agreement  is  obtained  with  the  previous  temperatures,  which  adds  further  support  to  rotational 
equilibriumL 


4.3.  Chemical  equilibrium 

The  question  of  chemical  equilibrium  was  addressed  using  the  kinetic  code  CHEMKIN  (Kee  et 
al.,  1986)  with  1 1  air  species  and  with  the  33  reactions  for  which  rates  were  selected  by  Park 
(1989).  In  this  code,  the  reverse  reaction  rates  are  computed  from  the  forward  rates  and  from 
thermodynamic  constants.  Since  the  CHEMKIN  thermodynamic  data  base  is  limited  to 
temperatures  below  5,000  K,  we  extended  the  range  up  to  10,000  K  using  the  thermodynamic 
functions  of  Liu  (1989, 1990). 

As  mentioned  in  Section  2.3,  the  maximum  temperature  reached  by  the  plasma  in  the 
induction  region  is  estimated  to  be  around  9,000  K.  At  the  point  where  the  spectral  measure¬ 
ments  are  made,  the  maximum  temperature  is  about  7,500  K  as  follows  from  the  measurements 
presented  in  Section  4.2.  The  flow  time  between  these  two  points  is  on  the  order  of  10  milli¬ 
seconds.  If  we  assume  that  the  plasma  in  the  coil  region  is  in  chemical  equilibrium  at  a 
temperature  of  9,000  K  and  undergoes  instantaneous  cooling  down  to  6,000  K,  then  the  time 
required  to  return  to  equilibrium  at  6,000  K  would  be  on  the  order  of  one  millisecond,  as  can  be 
seen  in  Fig  4.18.  Since  this  time  is  much  smaller  in  comparison  with  the  10  ms  plasma  char¬ 
acteristic  flow  time,  flow-induced  chemical  non-equilibrium  is  unlikely  in  this  plasma. 

Radial  diffusion  effects  and  rapid  cooling  can  induce  elevated  electron  number  densities  at  the 
plasma  edges.  However,  the  temperature  in  those  outer  regions  drops  below  5,000  K,  and 
therefore  contributes  only  a  negligible  percentage  to  the  radiative  power  emitted  along  a  plasma 
diameter. 
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5.  Spectral  emission  of  air  plasmas 


5.1.  Experimental  spectrum 

The  emission  spectrum  was  recorded  between  2,000  and  8,000  A  along  a  diameter  (5  cm)  of 
the  plasma.  This  wavelength  range  was  limited  by  the  poor  efficiency  of  the  detection  system  above 
8,000  A,  and  by  the  strong  absorption  of  oxygen  in  room  air  below  2,000  A.  Since  it  is  recorded 
along  the  plasma  diameter,  this  spectrum  includes  the  contribution  of  regions  at  temperatures  ranging 
from  300  to  7,500  K.  The  monochromator  entrance  and  exit  slits  were  respectively  set  at  60  and 
260  [im,  corresponding  to  a  trapezoidal  spectral  bandpass  of  base  3.5  and  top  2.2  angstroms.  The 
monochromator  scanning  speed  was  2.3  A/sec.  Data  were  acquired  using  the  lock-in  amplifier  at  a 
rate  of  5  Hz  with  an  averaging  time  constant  of  0.1  second. 

The  relative  spectral  response  of  the  detection  system  was  measured  with  a  deuterium  lamp 
(Optroiuc  Model  UV-40IR)  between  2,000  and  4,000  A,  and  a  quartz  halogen  tungsten  lamp 
(Optronic  Model  55A)  between  3,000  and  8,000  A.  Excellent  agreement  was  obtained  between  the 
two  calibrations  in  the  overlapping  region  3,000-4,000  A.  Absolute  intensity  calibration  was  ob¬ 
tained  with  the  quartz-halogen  tungsten  lamp  which  is  a  spectral  radiance  source  traceable  to  NIST 
standards.  With  these  two  high  quality  lamps,  calibration  uncertainties  were  limited  to  less  than 
10%,  with  the  maximum  in  the  extreme  regions  (2,000-3,000  and  7,000-8,000  A)  of  the  spectrum. 

A  high-bandpass  optical  filter  (cut-off  at  4,000  A)  was  used  when  recording  the  spectrum 
above  6,  000  A,  in  order  to  reject  the  second  order  structures  due  to  intense  near-UV  bands  such  as 
the  (0,0)  band  of  Nj  first  negative  and  the  (0,0)  to  (4,4)  bands  of  CN  violet  which  would  otherwise 
appear  strongly  in  the  vicinity  of  the  oxygen  triplet  at  7,774  A. 

The  experimental  spectrum  is  presented  in  Hg.  5.1.  As  mentioned  earlier,  it  includes  strong 
bands  of  NO,  Nj,  N2,  CN  and  lines  of  O,  N,  and  C.  The  spectral  region  5,000-7,000  A  is  not 
plotted  because  it  does  not  include  prominent  spectral  features  under  the  present  experimental 
conditions. 

5.2.  Computed  spectrum 

The  radiative  emission  spectra  presented  below  were  predicted  using  the  NEQAIR  code  (Park, 
1985a,  b).  The  required  input  data  are  the  chemical  composition  of  the  plasma,  the  kinetic  tempera¬ 
ture  of  the  heavy  species  (T)  and  of  the  free-electrons  (T^),  and  the  rotational  and  vibrational  tem¬ 
peratures  (Tp  T„)  of  the  molecular  species. 
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The  code  is  divided  into  two  parts,  an  excitation  and  a  radiation  module.  The  excitation 
module  predicts  the  population  distribution  of  the  electronic  states  in  the  Quasi-Steady-State  (QSS) 
approximation  using  a  collisional-radiative  model  for  the  populations  of  bound  electronic  states.  The 
radiative  module  determines  the  emission  (and  absorption)  by  performing  a  line-by-line  calculation  of 
the  radiative  contribution  of  the  electronic  systems  (Whiting,  1969). 

In  the  case  of  molecules,  the  contribution  at  a  particular  wavelength  X  of  a  rotational  transition 
between  two  vibrational  levels  v',  v"  belonging  to  different  electronic  levels  e’,  e"  is: 


/(X)= 


(5.1) 


ExaUdan 
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where  the  population  of  the  emitting  level,  the  probability  of  a  radiative  transition  from 

{e\y')  to  (c",v"),  Sjn  (Honi-London  factor),  the  probability  of  a  transition  from  J'  to  7",  AEy/j-  the 
difference  of  energy  between  the  two  levels,  and  d)(X)  the  lineshape  function  (Doppler,  coUisional, 
and  Stark  broadening). 

The  goal  of  the  present  work  is  to  test  the  radiative  module  of  NEQAIR  independently  of  the 
excitation  module.  This  is  possible  since  the  plasma  was  found  to  be  in  LTE.  The  QSS  computa¬ 
tions  of  the  excitation  module  were  therefore  by-passed,  and  the  population  of  the  electronic  levels 
replaced  by  a  Boltzmann  distribution: 
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=  (2/+l)exppf-^Mj.  (5.5) 

are  respectively  the  electronic,  vibrational,  and  rotational  partition  functions.  It  should  be  noted, 
however,  that  the  populations  of  most  radiating  states  are  predicted  to  be  within  10%  of  their  LTE 
value  when  using  the  QSS  approximation.  The  only  exception  is  the  upper  state  of  the  First  positive 
system  of  N^,  but  this  system  is  too  weak  to  be  measured  under  our  experimental  conditions. 

To  model  the  experiment,  the  plasma  was  divided  into  25  equally  spaced  concentric  rings,  and 
thermodynamic  properties  were  assumed  uniform  in  each  ring.  All  temperatures  (translational,  ro- 
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rational,  vibrational  and  electronic)  were  set  equal  to  the  LTE  temperature  measured  from  the 
7,774  A  oxygen  triplet.  In  each  ring,  the  chemical  equilibrium  composition  at  the  corresponding 
temperature  was  calculated  with  the  STANJAN  code  (Reynolds)  from  a  baseline  air  composition  of 
78.084%  N2,  20.946%  Oj,  0.934%  Ar,  and  0.033%  CO,.  This  small  amount  of  atmospheric  carbon 
dioxide  is  apparently  responsible  for  the  formation  of  CN,  whose  Violet  system  is  a  very  intense 
radiator,  and  of  C,  which  produces  a  strong  line  at  2,478.6  A. 

NEQAIR  was  used  to  predict  the  emission  spectrum  in  each  ring  between  2,0(X)  and  8,000  A, 
with  a  resolution  of  1,000  points  per  angstrom.  This  high  resolution  is  not  required  for  optically  thin 
computations,  but  is  necessary  for  absorption  calculations  in  order  to  resolve  the  rotational  lines, 
their  width  being  typically  0.01  A  under  the  conditions  of  this  experiment  The  radiative  transport 
equation  is  then  solved,  and  after  the  last  slab,  the  resulting  spectrum  is  convolved  with  the  experi¬ 
mental  slit  function  using  a  Fast-Fourier  transform. 

In  the  following  sections,  we  will  first  present  a  comparison  between  the  experimental  spec¬ 
trum  and  the  computations  performed  with  the  original  version  of  the  code,  NEQAIR- 1,  then  de¬ 
scribe  the  improvements  that  were  warranted  from  this  analysis,  and  finally  assess  the  quality  of  the 
revised  code,  NEQAIR-2. 

53,  Comparison  with  NEQAIR^l. 

The  comparison  with  the  predictions  of  NEQAJR-1  (see  Fig.  5.2)  revealed  that  several  features 
of  the  code  should  be  improved,  particularly  the  intensity  and  position  of  the  molecular  bands,  and 
that  additional  bands  should  be  considered.  As  may  be  seen  from  Figs.  5.1  and  5.2; 

(a)  In  the  shon  wavelength  range,  where  radiation  is  essentially  due  to  NO  bands,  the  numerical 
spectrum  fails  to  reproduce  both  the  intensity  and  the  shape  of  the  spectral  features. 

(b)  Between  2,950  and  2,975  A,  the  (4,2),  (3,1),  and  (2,0)  bands  of  Nj  (second  positive  system) 
arc  more  widely  spaced  than  in  the  experiment 

(c)  Between  3,000  and  4,000  A,  where  most  radiation  is  due  to  N,  (second  positive)  and  Nj  (first 
negative);  the  intensity  of  the  N,  bands  appears  to  be  a  factor  2  too  high  relative  to  the  Nt 
bands. 

(d)  The  intensity  of  the  violet  system  of  CN  (particularly  between  3,800  and  3,876  A)  is  strongly 
under  predicted  by  the  code. 

(e)  The  intensities  of  the  atomic  nitrogen  triplet  (7,423-7,442-7,468  A)  and  of  the  atomic  oxygen 
triplet  (7,773  A)  are  underestimated. 

This  initial  comparison  led  to  a  re-examination  of  the  radiative  model;  the  resulting  improve¬ 
ments  are  presented  in  the  next  section. 
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Figure  52.  Numerical  spectrum  {NEQAIR-l). 
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5.4.  Modeling  improvements. 


5.4.1.  Radiative  transition  probabilities: 

As  follows  from  Eq.  5.1,  the  emission  intensity  is  proportional  to  the  radiative  transition  prob¬ 
ability  This  quantity,  or  equivalently  the  square  of  the  electronic-vibrational  transition  mo¬ 

ment,  R^''  ,  was  approximated  in  NEQAlR-1  by  the  product  of  the  Franck-Condon  factor,  times 
the  square  of  the  electronic  transition  moment,  Since  the  variations  of  the  latter  with  in- 
temuclear  distance,  and  hence  with  the  vibrational  transition,  were  not  known  with  good  accuracy, 
Rg  was  taken  as  a  constant  in  NEQAIR-1  for  all  vibrational  bands  of  an  electronic  system.  However, 
recent  progress  in  computational  chemistry  has  resulted  in  accurate  predictions  of  the  electronic 
transition  moment  variations.  We  used  these  results  to  compile  reliable  transition  probabilities  for 
the  seven  electroiuc  band  systems  of  interest  in  our  temperature  and  spectral  ranges  that  were 
included  in  the  original  version  of  NEQAIR:  the  first  and  second  positive  systems  of  Nj,  the  first 
negative  system  of  Nj,  the  beta  and  gamma  systems  of  NO,  the  Schumann-Runge  system  of  Oj. 

A  detail  account  of  the  transition  probabilities  calculations  is  provided  below.  To  make  the 
results  applicable  to  a  broader  range  of  applications  including  optical  diagnostics,  quantitative 
spectroscopy  and  radiation  modeling,  the  computed  quantities  are  the  Einstein  A  coefficients,  the 
band  oscillator  strengths,  the  sums  of  the  squared  electronic-vibrational  transition  moments,  and  the 
Franck-Condon  factors. 

5.4.1.  a.  Definitions  and  previous  results 

The  key  parameter  from  which  transition  probabilities  are  calculated  is  the  square  of  the 
electronic-vibrational  transition  moment 

(r)R,(r)%..(r)dr)',  (5.6) 

where  v'  and  v"  are  the  vibrational  levels  in  the  upper  and  lower  electronic  levels  of  the  transition 
respectively,  'Fy'  and  'Fy"  are  the  corresponding  vibrational  wavefunctions,  r  is  the  intemuclear 
distance,  and  Re,  the  electronic  transition  moment  function  (ETMF),  is  an  average  of  the  electric 
dipole  moment  with  respect  to  the  electronic  wavefunctions. 

The  calculation  of  transition  probabilities  therefore  depends  on  the  determination  of  both  Re(r) 
and  the  vibrational  wavefunctions.  Although  the  latter  are  routinely  calculated,  the  electronic 
transition  moments  have  often  been  inferred  from  spectroscopic  intensity  measurements.  Such 
measurements  are  subject  to  difficulties  associated  with  gas  non-equilibrium,  presence  of  impurities, 
calibration,  and  low  signal  levels  for  weak  transitions.  Moreover,  the  most  Intense  bands  span  only  a 
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narrow  range  of  intemuclear  distance  r,  thus  experimental  band  strengths  give  very  httle  information 
about  the  variation  of  the  electronic  transition  moment  with  r.  Numerous  compilations  and  reviews 
of  transition  probabilities  based  on  these  measurements  are  nevertheless  available  (Nicholls,  1964  and 
1977,  Allen,  1966,  Schadee,  1967,  Kuz’menko  et  al,  1979). 

In  recent  years,  however,  electronic  transition  moment  functions  for  spin-allowed  transitions 
have  been  obtained  by  ab  initio  calculations  (Werner  et  al,  1984,  Rizzo  et  ai,  1988,  Langhoff  et  al, 
1988a,  1988b,  1991,  Bauschlicher  and  Langhoff,  1988).  Since  these  are  usually  in  very  good 
agreement  with  the  few  accurate  experimental  data  available,  it  is  then  reasonable  to  calculate  reliable 
transition  probabilities.  In  the  present  work,  we  combine  these  ab  initio  ETMFs  with  accurate 
spectroscopic  constants  to  produce  extensive  arrays  of  data  directly  usable  in  the  applications 
previously  mentioned. 


5.4.1. b.  Method  of  calculation 


The  r-centroid  approximation  (Fraser,  1954,  Nicholls  and  Jarmain,  1955)  has  been  widely  used 
to  calculate  the  integrals  of  Eq.  5.6.  Tlie  r-centroid,  defined  by 


rv-v" 


(5.7) 


has  such  properties  that,  under  a  wide  range  of  conditions,  the  following  approximation  holds  to  very 
good  accuracy  for  low  values  of  n  (Nicholls,  1974): 


-  f'F 

r"  =  ' 


■r''W-dr 


=  (rv'v")'' 


(5.8) 


Therefore,  when  the  ETMF  varies  with  the  intemuclear  distance  in  a  polynomial  fashion,  Eq.  5.6  may 
be  approximated  by 


=[R,(rvv")f  [jH',.(r)y,..(r)dr]',  (5.9) 

where  the  second  factor  on  the  RHS  is  the  square  of  the  vibrational  overlap  integral,  also  called  the 
Franck-Condon  factor.  The  importance  of  an  accurate  vibrational  analysis  is  clearly  seen  from 
Eq.  5.9,  since  Franck-Condon  factors  vary  by  several  orders  of  magnitude  from  band  to  band, 
whereas  electronic  transition  moments  usually  change  by  at  most  one  order  of  magnitude. 

The  r-centroid  approximation  gives  very  good  results  for  many  of  the  transitions  considered  in 
this  study,  with  the  possible  exception  of  NO  beta  for  which  the  ETMF  is  not  perfectly  fitted  by  a 
low  degree  polynomial.  The  method  also  requires  the  knowledge  of  r-centroids  and  Franck-Condon 
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factors,  which  have  not  always  been  calculated  firom  the  more  accurate  spectroscopic  data  now 
available. 

For  these  reasons,  we  have  instead  calculated  the  integrals  of  Eq.  5.6  directly.  This  required 
modification  of  an  existing  Franck-Condon  factor  program  (Zare,  1964).  Einstein  coefficients, 
absorption  band  oscillator  strengths,  and  sums  of  the  electronic-vibrational  transition  moments  are 
then  obtained  from  (Schadee,  1978); 


=  2.026  X  10^(v,.-)’  {*r'f .  (5.10) 


2-8„_. 


=  3.0376 X 10-* v^.-  ^  (fir" )’. 

She  (2-8,^.)  2-8„^. 


x(iir}  =(2-s,^.,j(2s+i)(Rr-) . 


(5.11) 

(5.12) 


provided  that  the  ETMFs  have  been  calculated  following  the  convention  of  Whiting  et  al  (1980).  In 
the  preceding  numerical  expressions,  Vv-v”  is  the  wavenumber  in  cm*l,  ''  the  electronic- vibrational 
transition  moment  in  atomic  units  (eao),  and  8i  j  the  Kronecker  delta,  which  is  equal  to  1  if  i=j  and  0 
otherwise;  (2S+1)  is  the  spin  multiplicity  and  A  the  projected  spin  angular  momentum  equal  to  0,  1, 
2, ...  for  L,  n.  A,...  states,  respectively. 

In  addition,  emission  band  oscillator  strengths  may  be  calculated  from  the  absorption  strengths: 


5.4.1. c.  Vibrational  analysis 


remis  _  2  ^q  a*  ^abs 
Jv'v’'  “  ^  c*  Jyv"' 
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(5.13) 


The  vibrational  wavefunctions  are  constructed  from  the  numerical  solution  of  the  radial 
Schrddinger  equation  of  nuclear  motion  (Cooley,  1961).  An  existing  program  (Zare,  1964)  was  used 
to  calculate  rotationless  Rydberg-Klein-Rees  potentials  based  on  Dunham  expansion  coefficients  T,y. 
The  program  requires  only  T/o  and  Yu  to  obtain 

G(v)  =  iyio(v+7y,  and  •  (5.14) 

i=l  1=0 

As  pointed  out  earlier,  special  attention  must  be  paid  to  the  spectroscopic  constants  Ty  that 
enter  the  construction  of  the  potentials,  Lavrov  (1991)  recently  published  a  review  of  such  constants 
for  O2,  Nt  and  NO.  For  N2,  however,  we  used  the  latest  results  of  Roux  and  Michaud  (1990)  and 
Roux  et  al  (1992).  For  N2''',  the  values  were  taken  from  Klynning  and  Pages  (1982),  and  for  CN, 
from  Rehfuss  et  al.  (1992)  and  Prasad  et  al.  (1992). 
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5.4.I.d.  Electronic  transition  moments 


N2  First  positive  system  B  —  A 

The  ETMF  from  Werner  et  al  (1984)  was  preferred  over  the  one  firom  Rizzo  et  al  (1988)  since 
it  yields  radiative  lifetimes  that  agree  better  with  the  experimental  values  summarized  in  Werner  et  al 
(1984),  Table  DC.  Piper  et  al  (1989)  inferred  Einstein  coefficients  by  cmnbining  their  experimental 
transition  moment  variation  with  the  radiative  lifetimes  measured  by  Eyier  and  Pipkin  (1983)  for  v'  = 
5  to  12;  this  procedure  enabled  them  to  determine  the  lifetimes  of  the  levels  v'  =  0  to  4,  and  these 
further  support  the  ETMF  from  Werner  et  al  (1984).  Fmally,  recent  lifetime  measurements  (1988) 
for  the  v'  =  0  level  yielded  a  value  of  13  ±  1  p.s  as  compared  to  13.5  p.s  when  using  the  ETMF  frcwn 
Werner  et  al  (1984)  vs.  9.8  jis  with  the  one  from  Rizzo  et  al  (1988).  The  quoted  uncertainty  on  the 
ETMF  is  5%,  and  thus  should  result  in  10%  accuracy  for  transition  probabilities. 

N2  Second  positive  system  C  —  B  -^11^ 

Werner  et  al  (1984)  calculated  an  ETMF  that  gives  excellent  agreement  with  the  measured 
radiative  lifetimes  presented  in  their  Table  VIII.  As  for  the  first  positive  system,  resulting  transition 
probabilities  should  be  accurate  within  10%. 

First  negative  system  B 

The  ETMF  calculated  by  Langhoff  and  Bauschlicher  (1988)  agrees  very  well  with  the 
experimental  Einstein  coefficients  of  Johnson  and  Fowler  (1970),  as  may  be  seen  in  Table  Vn  of 
Langhoff  and  Bauschlicher  (1988).  The  presently  calculated  transition  probabilities  are  in  even 
closer  agreement  because  updated  spectroscopic  constants  are  used  in  the  vibrational  analysis. 

NO  Beta  system  B  —  X^U.^ 

The  ETMF  of  the  NO  beta  system  is  more  difficult  to  determine  because  of  the  strong  coupling 
of  the  transition  upper  state  with  the  neighboring  C  state.  The  most  reliable  ETMF  was  obtained 
by  Langhoff  et  al  (1991).  The  resulting  lifetimes  are  somewhat  higher  than  the  experimental  values 
of  Gadd  and  Slanger  (1990)  for  v'  =  0  to  6.  However,  once  rescaled  by  the  ratio  of  the  experimental 
to  the  theoretical  lifetimes,  the  calculated  transition  probabilities  are  within  a  few  percent  of  the 
values  measured  by  Piper  et  al  (1991)  for  v'  =  0  to  5.  For  v'  =  6,  two  sets  of  experimental  transition 
probabilities  are  available  (Gadd  and  Slanger,  1990,  Piper  et  al,  1991).  While  the  two  sets  are  in 
good  agreement  (better  than  10%)  for  the  most  intense  bands,  larger  discrepancies  exist  for  the 
weaker  bands.  The  present  calculated  Einstein  coefficients,  once  rescaled  by  the  ratio  of  lifetimes, 
fall  however  between  the  two  sets.  Langhoff  et  al  (1991)  suggested  that  the  differences  may  be 
partially  attributed  to  calibration  problems  since  no  smooth  ETMF  can  match  both  the  lifetimes  and 
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the  transition  probabilities  of  Gadd  and  Slanger  (1990). 

For  levels  above  v'=6,  our  calculated  transition  probabilities  can  only  be  regarded  as 
approximate  because  of  the  increasing  influence  of  vibrational-rotational  interactions  (Langhoff  et  al, 
1991,  Gadd  etal,  1991). 

iVO  Gamma  system  A  —  X 

The  ETMF  used  here  is  the  (344)MRCI(7electrDns)  of  Langhoff  et  al  (1988).  The  calculated 
transition  probabilities  are  within  10%  of  those  measured  by  Piper  and  Cowles  (1986)  for  the  bands 
originating  from  v'  =  0, 1  and  2. 

O2  Schumann-Runge  system  B  —  X 

Transition  probabilities  for  the  O2  Schumann-Runge  are  obtained  with  the  semi-empirical 
ETMF  defined  by  Friedmann  (1990),  which  is  based  on  the  results  of  Allison  et  al  (1986).  The 
present  calculations  agree  very  well  with  the  experimental  mean  oscillator  strengths  of  Lewis  et  al 
(1986). 

CN  Violet  system  B  X 

The  transition  probabilities  of  the  CN  violet  system  are  computed  using  the  electronic  transition 
moment  of  Bauschlicher  and  Langhoff  (1988). 


5.4.1.e.  Results 

Tables  of  Einstein  A  coefficients,  absorption  band  oscillator  strengths,  and  sums  of  the  elec¬ 
tronic-vibrational  transition  moments  are  given  for  each  of  the  systems  previously  mentioned. 
Franck-Condon  factor  arrays  hare  also  been  computed  so  that  comparisons  can  be  made  with  other 
sets  of  experimental  transition  probabilities  obtained  with  different  Franck-Condon  factors.  Rnally, 
tables  of  wavelengths  and  r-centroids  have  also  been  prepared.  All  tables  are  available  in  the  com¬ 
puter  diskette  attached  to  this  report.. 

The  1,500  or  so  updated  transition  probabilities  that  were  introduced  into  the  new  version  of 
the  NEQAIR  code  have  resulted  in  intensity  changes  of  about  -10%  for  the  first  negative  system  of 
Nj  and  the  second  positive  systems  of  N2,-15%  for  the  first  positive  system  of  N^,  -25%  for  NO 
beta,  +300%  for  NO  gamma,  +250%  for  the  Schumann-Runge  system  of  Oj,  and  a  factor  10  in¬ 
crease  for  the  Violet  system  of  CN  (these  values  are  for  the  global  power  emitted  between  2,000  and 
20,000  A  by  a  plasma  in  LTE  at  8,000  K,  and  are  virtually  constant  over  at  least  the  range  4,000- 
8,000  K). 
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5.4.2.  Spectroscopic  constants: 


Predicting  line  positions  accurately  is  of  special  importance  when  trying  to  match  experimental 
measurements  to  gain  information  on  populations  or  temperatures,  or  when  absorption  calculations 
are  required.  Line  positions  are  determined  from  the  difference  of  energy  between  the  levels  of  the 
transition, 

Ar  =  Ar,  +  AG(v)+AF„(7)  (5.15) 

where  T^,  G(v),  and  FJJ)  are  the  electronic,  vibrational,  and  rotational  term  energies  respectively. 
G(v)  is  itself  defined  as 

G(v)  =  0)^  (v + 1)  -  {y+jf  +  ( V + +  a),z,  ( V  + 1)“  +. . .  (5. 16) 

and  FJJ)  as 

F,U)  =  B,J(J  +  \)-D,J^iJ  +  l)^+...,  (5.17) 

where 

=  fl*(v+y)-a,(v+i)  +Y^(v+^)  +5,(v  +  i)  +...,  (5.18) 

D.=D.{v+i)+P,{v+if+...  .  (5.19) 

Improved  electronic  term  energies  (T^)  and  Klein-Dunham  coefficients  (©^,  to^*,...,  a,,  D^, 

P,,...),  based  on  accurate  data  from  Fourier  transform  spectroscopy,  have  been  included  into 
NEQAlR-2;  they  are  tabulated  and  referenced  in  Appendix  Cl.  Since  these  spectroscopic  constants 
were  obtained  by  polynomial  fits  of  experimental  line  positions,  it  is  essential,  when  using  such  de¬ 
velopments,  to  include  the  same  number  of  polynomial  terms  as  in  the  original  fit  NEQAIR-1,  which 
was  limited  to  3  terms  in  each  development,  was  modified  to  include  all  necessary  coefficients.  The 
resulting  effect  is  quite  drastic  for  the  vibrational  bands  of  the  second  positive  system  of  N, 
(Fig.  5.3);  the  large  shift  of  about  30  A  in  band  positions  in  NEQAlR-1  was  due  to  the  neglecting  of 
(D^z,  in  Eq.  5.16. 
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Figure  S3  Effect  of  improved  spectroscopic  constants  on  the  N2  Second  positive  system. 


XlA] 


Figure  5  Ja.  Experimental  spectrum 


XlA) 


Figure  5  Jb.  Numerical  spectrum  (NEQAIR-I). 


X[A) 

Figure  S.3c.  Numerical  spectrum  (NEQAIR-2). 
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5.43.  Additional  bands: 


The  spectral  region  below  2,800  A,  where  the  electronic  systems  of  NO  are  the  dominant  ra¬ 
diators,  is  of  particular  interest  for  ultra-violet  signatures.  It  appeared  in  the  course  of  this  work  that 
the  numerical  spectrum  failed  to  reproduce  the  experimental  one,  even  after  updating  the  u-ansition 
probabilities  for  NO  Beta  and  Gamma  (see  Fig.  5.4a).  We  therefore  considered  additional  electronic 
systems  of  NO,  including  the  Delta,  Epsilon,  Beta  prime,  and  Gamma  prime.  The  positions  of  the 
emitting  levels  for  these  ttansitions  are  sketched  in  fig.  5.5.  The  electronic  level  C  (upper  of  the 
Delta  transition),  is  strongly  predissociated  above  v'=0,  J’>4.5  (see  for  example  Gadd  et  al.,  1991) 
due  to  interferences  with  the  a"*!!  state,  but  the  other  levels  are  not  subject  to  predissociation.  We 
present  first  the  addition  of  the  non-predissociated  levels,  and  then  the  contribution  of  NO  Delta. 

5.4.3. a.  NO  Epsilon  (D-^  X): 

The  radiative  transition  probabilities  for  this  system  were  obtained  from  the  electronic  transi¬ 
tion  moment  computed  by  De  Vivie  and  Peyerimhoff  (1988),  and  from  potential  curves  based  on  the 
spectroscopic  constants  of  Huber  and  Herzberg  (1979)  and  Amiot  (1982).  However,  while  compu¬ 
tational  chemistry  usually  predicts  the  correct  variation  of  the  electronic  transition  moment  with  in- 
temuclear  distance,  experimental  studies  produce  more  accurate  values  of  its  absolute  magnitude. 
For  this  reason,  the  radiative  lifetime  of  the  v’=0  state  was  matched  to  the  radiative  lifetime  measured 
by  Tsukiyama  et  al.  (1988),  and  the  transition  probabilities  were  normalized  accordingly. 

As  shown  in  Fig.  5.4b,  the  NO  Epsilon  bands  contribute  significantly  to  the  radiation  in  the 
range  2,000-2,500  A  .  In  particular,  the  “stair-case”  structure  of  the  experimental  spectrum  in  that 
region  is  numerically  recovered,  at  least  partially,  with  the  introduction  of  NO  Epsilon.  Although 
further  studies  are  needed  to  predict  band  intensities  with  better  accuracy,  NO  Epsilon  appears  as  a 
major  source  of  radiation  in  the  ultra-violet 

5.4.3. b.  NO  Beta  Prime  (B'— >  X)  and  Gamma  prime  >  X): 

The  radiative  teansition  probabilities  for  NO  Beta  prime  and  Gamma  prime  were  obtained  from 
the  electronic  transition  moments  of  De  Vivie  and  Peyerimhoff  (1988),  and  from  the  spectroscopic 
constants  of  Huber  and  Herzberg  (1979)  and  Amiot  (1982).  These  ffansition  probabilities  are  about 
one  order  of  magnitude  smaller  than  those  of  the  Delta  or  Epsilon  transition.  Moreover,  since  the 
emitting  levels  lie  about  1  eV  higher  in  energy  than  the  C  or  D  states,  their  population  is  much 
smaller  than  those  of  the  C  and  D  states,  and  the  bands  lie  farther  in  the  vacuum  ultra-violet.  For 
these  reasons,  the  NO  Beta  prime  and  Gamma  prime  bands  produce  a  very  small  amount  of  radiation 
in  the  wavelength  and  temperature  ranges  considered  here.  The  missing  structures  in  the  computed 
spectrum  of  Fig.  5.4b  cannot  therefore  be  attributed  to  these  systems. 
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Figure  5.4.  Influence  of  additional  NO  bands. 


Fig  S.4a.  NO  Beta  and  Gamma.  Optically  thin 
computations. 


Fig  5.4d.  NO  Beta,  Gamma,  Delta,  Epsilon  with 
predissociation  of  C  state.  Optically  thin 
computations. 


Fig.  5.4b.  Addition  of  NO  Epsilon.  Optically  thin 
computations. 


Fig.  5.4e.  NO  Beta,  Gamma,  Delta  and  Epsilon 
with  predissociation  of  the  C  state. 
Including  self  absorption. 


Fig.  S.4c.  NO  Beta,  Gamma,  Epsilon  and  (bottom 
curve)  NO  Delta  without  predissociation. 
Optically  thin  computations. 
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Figure  5^.  NO  energy  level  diagram. 


5.43.C.  NO  Deltaic^  X): 

The  radiative  transition  probabilities  for  NO  Delta  were  determined  from  the  electronic  tran¬ 
sition  moments  of  Cooper  (1982)  and  from  the  spectroscopic  constants  of  Huber  and  Herzberg 
(1979)  and  Amiot  (1982). 

As  mentioned  earlier,  the  emitting  C-level  is  strongly  predissociated.  Indeed,  if  such  were  not 
the  case,  the  computed  intensity  of  the  Delta  band  would  be  much  too  strong,  as  may  be  seen  in 
Fig.  5.4c.  It  is  interesting  to  note  however  that  the  position  of  the  missing  features  of  the  numerical 
spectram  coincide  with  the  strongest  Delta  bands. 

Several  authors  have  measured  (Tsuidyama  et  al.,  1987  and  1988,  Rottke  and  Zacharias,  1985, 
Benoist  dAzy  et  al.,  1975,  Hildda  et  al.,  1987,  Hart  and  Bourne,  1989)  or  computed  (DeVivie  and 
Peyerimhoff,  1990)  the  predissociative  lifetimes  of  individual  rotational  levels  of  the  C  state,  and 
have  obtained  values  on  the  order  of  30  ns  for  the  non-predissociated  levels  (v=0,  r<5),  and  between 
10  ns  and  0.03  ns  for  the  predissociated  ones.  On  the  other  hand,  the  kinetic  collision  time  (i.e.  the 
inverse  of  the  collision  frequency)  by  heavy  particle  impact  is  on  the  order  of  0.3  ns  at  7,000  K  and 
1  atm.  We  argue  below  that,  in  our  high  pressure  (1  atm)  and  high  temperature  experimental 
conditions,  collisional  processes  compete  with  predissociation  to  maintain  some  population  in  the 
predissociated  levels. 

A  complete  coUisional-radiative  model  would  be  required  to  predict  the  population  of  all  rota- 
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tional  levels  in  the  C  state.  Ideally,  the  following  mechanisms  would  be  included: 

(a)  Rotational  relaxation  within  C  state  by  heavy  particle  impact  (rotational  transitions  occur 
through  momentum  transfer,  and  electrons  are  inefficient  in  transferring  momentum^^) 

(b)  Collisional  excitation  and  quenching  of  the  C  state  by  heavy  panicles  and  electrons. 

(c)  Recombination  of  N  and  O  atoms,  or  NO  and  electrons  to  the  C  state. 

(d)  Radiative  cascading. 

(e)  Predissociation  of  the  C  state. 

Such  a  naodel  is  unfortunately  out  of  reach  presently  for  lack  of  knowledge  of  most  cross  sec¬ 
tions  for  the  above  processes.  However,  Callear  and  Pillin  (1970)  and  Lahmani  et  al  (1981)  ob¬ 
served  fast  collisional  exchange  between  the  C  and  the  D  state  in  the  presence  of  argon.  This  is  not 
surprising  in  light  of  the  proximity  of  the  two  levels  (the  difference  of  energy  between  C,  v=0,  J=0 
and  D,  v=0,  J=0  is  about  900  cm-‘,  or  about  a  tenth  of  an  eV).  In  both  studies,  the  collision  rate  in 
the  presence  of  argon  was  found  to  be  close  to  the  kinetic  collision  rate.  Le  Duff  et  al.  (1986)  ob¬ 
served  the  same  phenomenon  in  the  presence  of  other  atomic  gases  (rare  gases).  On  the  other  hand, 
molecular  gases  were  found  to  be  much  less  efficient  quenchers  (Callear  and  Pilling,  1970).  In  our 
experiment,  however,  the  plasma  in  the  region  of  intense  radiative  emission  is  mainly  composed  of 
atomic  species  (N  and  O)  — the  total  mole  fraction  of  atoms  is  over  0.90  at  7,000  K. 

Therefore,  considering  only  predissociation  and  collisions  with  the  D  state  by  heavy  particle 
impact,  i.e. 

NO(C,  V,  J)+  M<-^NO(D,v’,r)  +  M, 

where  M  is  a  heavy  particle  (atom),  the  population  evolution  of  a  given  rotational  level  J  in  the  C 
state  is  given  by: 

(5.20) 

where  is  the  collisional  excitation  rate  of  the  J  level,  its  reverse,  and  the  concentration 
of  heavy  particles. 

Using  detailed  balance,  the  reverse  rate  can  be  expressed  as 

{n 

(5.21) 

Assuming  that  the  D  level  which  is  not  subject  to  predissociadon  (see  for  instance  Callear  and  Pilling, 
1970)  remains  in  LTE,  and  assuming  quasi-steady  state,  Eq.  5.20  becomes: 


(5.22) 
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or,  assuming  unit  efficiency  for  heavy  particle  collisions  as  suggested  by  Callear  et  Pilling  (1970): 


The  collision  frequency  was  approximated  by  (Vincenti  and  Kruger,  1965): 


colluions 
^  '^colUsiotis 


(5.23) 


=  10'®  P[atm], 


(5.24) 


and  the  values  kprtd  were  obtained  from  a  fit  of  the  few  available  experimental  and  theoretical  pre¬ 
dissociation  data  (Hg.  5.6).  Because  of  spin  interaction,  each  rotational  level  J  is  split  into  Q=l/2 
and  Q=3/2  substates.  Fig  5.6  shows  that  the  £2=1/2  components  have  a  much  shorter  predissociative 
lifetime  than  the  £2=3/2  component  above  the  predissociation  limit  For  this  reason,  all  emission 
from  £2=1/2  above  J=4.5  was  neglected. 
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Figure  5.6.  Predissociative  lifetimes  of  the  rotational  levels  of  the  C  state  of  NO  (v'=0). 


The  numerical  spectrum  predicted  by  NEQAIR-2,  with  this  simplified  coilision-predissociation 
model  included,  is  presented  in  Rg  5.4d.  All  previously  missing  features  in  the  computed  spectrum 
are  numerically  recovered  with  the  simple  model  developed  above.  The  absolute  magnitude  of  the 
new  spectrum  is  however  a  little  larger  than  in  the  experiment,  but  self-absorption  by  the  plasma 
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should  be  considered.  As  may  be  seen  in  Fig.  5.4e,  including  self-absorption  leads  to  excellent 
agreement  in  shape  and  absolute  intensity  between  the  experimental  and  numerical  spectra. 

Of  course,  more  effort  is  needed  to  predict  accurately  the  radiation  due  to  NO  Delta  and 
Epsilon.  Among  the  most  important  data  required  for  the  collision-predissociation  model  are  better 
estimates  of 

(a)  predissociation  rates  for  the  rotational  levels  of  the  C  state, 

(b)  collisional  excitation  cross-sections  between  the  C  and  the  D  states, 

(c)  radiative  transition  probabilities  for  the  Delta  and  Epsilon  transitions. 

5,5.  Comparison  with  NEQAIR-2. 


The  final  comparison  between  experiment  and  improved  computations  is  shown  in  Fig.  5.7. 
Much  better  agreement  is  now  obtained  for  the  shape,  position,  and  intensity  of  the  molecular  bands. 
The  additional  Delta  and  Epsilon  bands  of  NO  have  significant  effect  in  the  ultra-violet  region  which 
is  now  well  reproduced  numerically.  The  addition  of  NO  Beta  prime  and  Gamma  prime  bands  has  a 
negligible  incidence  on  the  spectrum  presented  in  Fig.  5.7,  but  would  have  more  importance  in  the 
spectral  region  below  2,000  A,  where  the  most  intense  bands  lie,  or  at  higher  temperatures  where  the 
emitting  levels  are  more  populated.  The  position  of  the  N2  Second  positive  bands  has  been 
improved;  their  intensity,  too  high  by  a  factor  2  in  NEQAIR-1,  was  corrected  by  the  addition  of  the 
homonuclear  factor  in  partition  functions  calculations  as  described  by  Moreau  et  al  (1992).  The  CN 
Violet  bands  now  appear  with  their  correct  intensity.  Finally,  the  intensities  of  the  nitrogen  and 
oxygen  triplets  are  closer  to  the  measurements.  This  change  is  not  due  to  any  physical  improvement, 
but  rather  to  a  higher  numerical  resolution;  NEQAIR-1  used  constant  wavenumber  spacing  instead  of 
constant  wavelength  spacing  and  therefore  few  spectral  points  were  considered  at  long  wavelengths. 

Some  spectral  features  still  appear  underpredicted  by  NEQAIR-2  (Nj  (0,0)  band  at  3,911  A, 
some  Nj  bands,  and  the  O  and  N  triplets).  These  differences  (about  20  %)  are  however  within 
modeling  and  experimental  uncertainties.  In  any  case,  predicting  the  intensity  of  the  oxygen  triplet  at 
7,773  A  with  10  %  accuracy  would  require  the  temperature  to  be  measured  with  0.5  %  accuracy. 
The  remaining  intensity  differences  are  therefore  not  significant 
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Figure  5.7b.  Numerical  spectrum  {NEQAIR-2). 


6.  Radiative  source  strength  of  air  plasmas 


For  atmospheric  pressure  air  at  temperatures  above  9,000  K,  measurements  of  the  radiation 
source  strength,  or  local  emissivity,  have  been  reported  by  several  authors,  including  Page  and 
Arnold  (1964),  Nerem  (1964),  Schiciber  et  al.  (1973),  Ogurtsova  et  al.  (1968),  and  Devoto  et  al. 
(1978).  These  experiments  were  conducted  in  shock  mbes  (Page  et  al.,  Nerem),  or  in  arc  dischai^es 
(Schreiber  et  al.,  Ogurtsova  et  al.,  Devoto  et  al.).  The  emitted  radiation  was  recorded  by  thermopiles 
and  calibrated  with  tungsten  lamps.  Temperatures  were  determined  under  the  LTE  assumption  firom 
the  absolute  intensity  of  atomic  oxygen  and  nitrogen  lines. 

In  addition,  several  theoretical  nodels  have  been  developed  over  the  past  decades  (Kivel  and 
Bailey,  1957,  Meyerott  et  al.,  1960,  Breene  et  al.,  1963,  and  Park,  1985).  A  comparison  between 
the  measurements  and  the  predictions  from  of  the  codes  is  presented  in  Fig.  6.1.  The  curve  labeled 
as  NEQAIR-2  includes  the  changes  described  in  the  previous  sections. 


Figure  6.1  Volumetric  radiative  source  strength  of  LTE  air  at  atmospheric  pressure:  previous  results  and  en¬ 
hanced  NEQAIR-2  computations. 


The  measurements  conducted  during  the  course  of  this  work  cover  the  range  4,000  to  8.000  K 
and  are,  to  our  knowledge,  the  first  measurements  in  that  temperature  range.  The  radiation  source 
strength,  or  local  emissivity,  as  a  function  of  lateral  position  was  measured  as  described  in  Section  3. 
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The  lateral  emission  profile  was  then  corrected  for  plasma  self-absorption  using  the  following 
procedure.  The  plasma  was  modeled  using  the  LTE  temperature  profile,  corrected  for  self 
absorption,  measured  in  Section  4.  The  lateral  spectrally  integrated  emission  profile  was  computed 
for  both  optically  thin  and  self-absorbed  cases  using  the  enhanced  version  of  the  NEQAIR  code, 
NEQAlR-2.  The  measured  profile  of  spectrally  integrated  emission,  thereby  corrected  for  self 
absorption  effects,  was  then  Abel-inverted  and  reduced  to  volumenic  radiative  source  strength.  It 
should  be  noted  here  that  the  self-absorption  correction  is  on  the  order  of  about  10%.  Finally,  the 
spectrally  integrated  radiation  was  plotted  against  the  self-absorption  corrected  LTE  temperature 
from  the  7,773  A  oxygen  line.  The  resulting  radiation  source  strength  profile  is  presented  in  Figure 
6.2,  where  it  is  compared  to  the  predictions  of  the  aforementioned  codes.  In  all  cases,  the  spectral 
range  extends  from  2,000  to  20,000  A. 

While  the  early  models  overestimate  the  radiative  source  strength,  the  enhanced  code 
NEQAIR-2  agrees  closely  with  the  present  measurements. 


Temperature  [K] 


Figure  6.2  Volumetric  radiative  source  strength  of  LTE  air  at  atmospheric  pressure:  comparison  between  our 
experimental  results  and  the  predictions  of  numerical  codes. 
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7.  Conclusions 


Optical  diagnostics  were  conducted  in  an  atmospheric  pressure  plasma  torch  operating  on  air  in 
the  purpose  of  studying  fundamental  plasma  processes.  The  first  task  consisted  in  determining  the 
thermodynamic  state  of  the  air  plasma  through  atomic  and  molecular  emission  spectroscopy.  The 
plasma  was  found  to  be  close  to  Local  Thermodynamic  Equilibrium,  and  composition  and 
temperature  profiles  were  determined.  The  spectrum  emitted  by  the  air  plasma  over  the  spectral 
range  2,000  to  8,000  A  was  then  compared  to  the  LTE  predictions  of  the  NEQAIR  code.  This 
comparison  led  to  a  number  of  improvements  to  the  radiative  model  of  the  code.  In  particular,  the 
transition  probabilities  and  spectroscopic  constants  of  the  main  radiating  systems  were  updated,  and 
additional  electronic  systems  were  included  in  the  code.  A  data  base  was  constituted  for  several 
band  systems  including  the  First  and  Second  positive  systems  of  ,  the  First  negative  system  of 
Nt,  the  Beta,  Gamma,  Delta,  Epsilon,  Beta  prime  and  Gamma  prime  systems  of  NO,  the 
Schumann-Runge  system  of  O2,  and  the  Violet  system  of  CN.  Among  the  new  electronic  systems 
introduced  in  the  code,  NO  Delta  and  Epsilon  were  found  to  be  especially  strong  radiators  in  the 
ultra-violet.  A  collision-predissociation  model  for  the  NO  C-state  from  which  the  Delta  transition 
originates  is  however  necessary,  and  a  simplified  such  model  has  been  presented.  After  the  owdeling 
changes,  excellent  agreement  in  shape  and  absolute  intensity  is  obtained  between  the  expetimental 
and  numerical  spectra.  The  radiative  module  of  NEQAlR-2  now  predicts  radiative  intensities  within 
20%  accuracy  for  most  bands  over  at  least  the  spectral  range  2,000-8,000  A.  Fmally,  the  radiation 
source  strength,  or  local  emissivity,  was  measured  for  the  first  time  over  the  temperature  range  4,000 
to  8,000  K.  The  predictions  of  the  enhanced  NEQAIR  code  are  found  to  be  within  modeling  and 
experimental  error  of  our  source  strength  measurements. 
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